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1 .   INTRODUCTION 

Many  approaches  to  the  computer  aided  design  problem  are 

being  studied.   To  date  these  efforts  have  not  resulted  in  the  general 

2 
availability  of  automatic  circuit  design  programs  .  If,  however,  one 

subdivides  the  electronic  circuit  design  process  somewhat,  one  of  the 
principle  facets  which  emerges  is  the  need  for  a  through  analysis 
capability.  With  such  a  tool  available  to  the  designer,  electronic 
system  response  may  be  simulated  over  the  entire  space  of  component 
variations  and  operating  conditions. 

One  of  the  fundamental  steps  in  designing  electronic  cir- 
cuits is  the  choice  of  the  "static  d.c."  or  quiescent  operating  point. 
Its  optimization  is  often  an  iterative  process.  An  example,  from 

logic  circuit  design,  occurs  when  designing  emitter  coupled  logic 

3 
circuits  such  as  those  of  Figure  1.   The  typical  design  process 

requires  that  a  number  of  interdependencies  be  met  in  order  to  obtain 
the  nonsaturating  output  characteristics  of  Figure  2.   The  process 
begins  with  the  selection  of  the  voltage  levels  representing  each 
binary  state.  The  on  and  off  state  base-to-emitter  voltages  may  be 
determined  by  superimposing  the  source  line  upon  a  set  of  the  tran- 
sistor 's  input  characteristic  curves.   The  output  levels  are  next 
determined  for  various  loadings,  such  as  the  next  stage  of  Figure  1A, 
or  the  output  level  restoring  emitter  followers  of  Figure  IB,  using 
the  transistor's  output  characteristic  curve  family.   This  procedure 

is  then  repeated  for  varying  input  and  output  conditions  to  insure 

5 
the  convergence  of  the  logic  levels  to  the  desired  voltage  levels. 
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Figure  1.  Two  Forms  of  Emitter  Coupled  Logic  (ECL) 
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The  amount  of  calculation  required  for  effective  design  of  an  emitter 
coupled  logic  block  can  "be  substantial  due  to  the  numerous  iterations 
which  may  be  needed  to  prove  satisfactory  operation  under  varying 
tolerance  conditions.   The  point  being  made  here  is  that  the  circuit 
designer  could  save  a  significant  amount  of  computation  time  by 
employing  a  flexible  computer  program  whose  output  could  be,  for 
example,  the  emitter  coupled  logic  circuit's  output  voltage  transfer 
characteristic . 

The  above  "static"  design  discussion  employed  graphical 
analysis  which  has  been  a  standard  circuit  design  tool  for  many 
years.   Introductory  electronic  circuit  texts  introduce  this  approach 
as  being  the  best  way  to  handle  nonlinear  circuit  elements  and  indeed 
it  is  for  simple  exploratory  analysis.  However,  for  computer-aided 
design,  such  as  is  possible  with  a  parameter  variation  program,  it  is 
necessary  to  enter  the  transistor's  characteristics  in  the  form  of 
equations.   These  equations  are  known  as  the  transistor  modeling 
equations  and  may  take  many  forms. 

Some  of  the  models  currently  used  are:  Ebers-Moll, 
Piecewise  Linear,  Admittance  or  Y  Parameter,  Hybrid  or  H  Parameter, 
Z  matrix,  and  the  M  Parameter  Model.   Each  model  evolved  for  the 
purpose  of  simulating  some  portion  of  a  transistor's  range  of  opera- 
tion. No  analysis  or  simulation  program  can  be  effective  unless  the 
transistor  model  used  accurately  represents  the  transistor's  voltage- 
current  behavior  throughout  the  region  of  operation  in  the  simulated 
circuit.   For  this  reason,  the  basis  of  each  of  the  above  models  will 


be  discussed  in  section  2  in  order  that  they  may  be  compared  with  a 

7 
new  model  related  to  the  Z  matrix  model  but  including  the  additional 

constraints : 


BETA  =  A  +  A  (In  ^  )  +  A   (in  ■=$■  f 
U  0  0 


VBE  =  B0  +  Al  (in  TQ   >  +  A2  <ta  |  )2 


These  additional  constraints  are  formed  by  the  Linear 

Regression  program  of  Appendix  B  whose  input  data  points  are  taken 

from  measurements  or  from  manufacturer's  data  sheets.   This  Regression 

Model  thereby  enables  the  designer  to  employ  linear  equation  solution 

methods  to  obtain  a  nonlinear  solution  to  his  circuit  problem.   Thus 

the  quantities  BETA  and  V^  are  no  longer  restricted  to  constant  values 

Bill 

as  is  the  usual  case. 

The  usefulness  of  the  Regression  Model  rests  upon  the  ease 
of  its  application  and  upon  the  accuracy  of  its  predictive  capability. 
In  an  effort  to  satisfy  the  first  goal,  the  Regression  Model  sub- 

o 

routine  written  by  Guth  was  converted  from  IBM  1620  compatible 
language  to  FORTRAN  IV  for  use  on  the  IBM  360  Model  75  as  a  part  of 

this  thesis  project. 

9 
In  the  process  of  conversion  a  new  version  of  ECAP  became 

available  which  provided  true  worst  case  solutions.   This  new  version 
of  ECAP  was  modified  to  work  with  an  updated  version  of  Guth's  sub- 
routine, a  listing  of  which  is  given  in  Appendix  B  of  this  thesis. 


The  complete  program  is  now  called  NLAP  (Nonlinear  Analysis  Program) 
and  is  contained  in  the  program  library  at  the  Department  of  Computer 
Science,  University  of  Illinois.  The  KLAP  calling  sequence  is  given 
in  Appendix  E,  together  with  the  input  and  output  of  a  test  program. 

The  final  goal  posed  for  this  thesis  was  to  evaluate  the 
possibility  of  extending  the  useful  range  of  the  Regression  Model 
from  the  normal  active  region  to  include  operation  in  the  cutoff 
region  as  is  encountered  in  emitter  coupled  logic  circuit  operation. 
The  evaluation  of  this  extension,  together  with  a  discussion  of  the 
IBM  36O  ECAP  network  solution  method,  will  be  presented  in  section  3« 

The  application  of  NLAP  to  the  solution  of  transistor  logic 
circuits  is  discussed  in  section  k.      The  required  topological  codes 
are  given  for  the  various  types  of  branch  level  statements.  Two 
illustrative  examples  are  presented  as  a  further  guide.   It  is 
recommended  that  the  user  consult  sections  3«^+  through  3*6  before 
undertaking  use  of  the  NL  solution  option  discussed  in  section  h.3- 


2.   THE  DEVELOPMENT  OF  THE  REGRESSION  MODEL  AND 
ITS  RELATIONSHIP  TO  CONVENTIONAL  CIRCUIT  MODELING 


2.1  Basic  Modeling  Elements 

The  basic  components  for  an  electronic  model  are  the  ideal- 
ized pure  resistance  (R),  pure  inductance  (L),  and  pure  capacitance  (C) 
elements  of  linear  and  reciprocal  network  theory.   Their  elemental 
voltage -current  relationships  are  simply  Ohm's  law.   To  form  a  compact 
set  of  elements,  it  is  necessary  to  add  the  controlled  source  for 
modeling  of  the  amplification  process  and  the  ideal  junction  diode  for 
simulating  variable  amplitude,  unilateral  transmission. 

Dependent  current  and  voltage  controlled  sources  are  of 
four   possible  types;  voltage  controlled  voltage,  voltage  controlled 
current,  current  controlled  voltage,  and  current  controlled  current. 
The  internal  impedances  of  these  sources  are  identical  with  their 
independent  counterparts  of  linear  network  theory.   The  direction 
of  output  current  flow,  or  of  output  voltage  polarity,  is  dependent 
upon  the  driving  source  and  may  be  of  either  polarity.   The  two  port 
volt age -current  relationships  of  each  controlled  source  type,  together 
with  their  schematics,  are  summarized  by  Figure  3« 

The  ideal  junction  diode  has  the  voltage -current  relationship: 

.   /  qV/MKT  .  v 
i  =  iQ  (eH  '    -1). 

This  equation  derives  from  a  solution  of  the  one -dimensional  diffusion 

11 
equation.    The  factor  M  in  the  denominator  of  the  exponent  facili- 
tates adjusting  the  theoretical  V-I  relationship  to  match  that  observed 
physically.  Use  of  this  ideal  diode  equation,  together  with  the  other 
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(a)  INPUT  CURRENT  CONTROLS  OUTPUT  CURRENT 
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(b)  INPUT  VOLTAGE  CONTROLS  OUTPUT  CURRENT 


(c) INPUT  CURRENT  CONTROLS  OUTPUT  VOLTAGE 


(d)  INPUT  VOLTAGE  CONTROL  OUTPUT  VOLTAGE 
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Figure  3*  Dependent  Sources  and  Their  Voltage -Current  Characteristics 


four  basic  modeling  elements,  plus  an  independent  current  or  voltage 

12 
source  element   ,  are  sufficient  for  most  modeling  work  whether  the 

circuit  model  be  derived  from  basic  physics  or  from  the  circuit 

"black  box"  approach. 


2.2  Small  Signal  Parameter  Models 

A  general  network  or  "black  box"  may  be  described  in  terms 
of  its  two  port  parameters.  At  each  port  the  voltage  and  input  current 

are  measured,  providing  four  variables  for  network  description.  Any 

13 
two  of  these  quantities  form  a  basis  "  for  the  solution  space  of  the 

network  and,  therefore,  they  completely  characterize  the  network's 

behavior.   The  two  port  parameters  may  be  established  from  a  graphical 

plot  by  following  the  definitions  established  during  their  derivation. 

One  pair  of  general  relationships  between  the  two  port 

voltages  and  current  is: 


ei  =  t{ii>  V 


e2  =  g  (!-(_,  i2) 


These  general  functional  relationships  will  in  general  be 
nonlinear.   This  creates  the  need  for  their  transformation  into  linear 
functions  for  use  in  linear  network  analysis.   Expanding  in  a  Taylor 
series,  with  two  independent  variables,  obtains: 
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el  (V  i2)  = 


>  (l_,  Ip)  +  °el 

^1 


(in  - 


i_  =  I. 


I,)  +  3el 
1    ^ 


(i2  -  I2) 


I«  =  I- 


+  higher  order  terms 


(V  12)  = 


.2(lr  I2)+H 


(in  " 


=  I 

: 
=  i. 


d2  - 12) 


i„  =  I, 


+  higher  order  terms 

provided  e,  (i.   ip)  and  ep  (i  ,  i  )  are  analytic  functions.   The 
crucial  point  is  that  there  are  two  ways  for  these  Taylor  expansions 
to  reduce  to  linear  functions.   First,  if  the  basic  functional  rela- 
tionships are  linear,  then  the  value  of  all  derivatives  of  higher 
order  than  one  are  zero.   Second,  if  the  operating  point  Q  (i  ,  i  ) 
is  selected  such  that  the  deviations  (i  -  I  )  and  (i  -  I  )  are  small 
(thereby  restricting  the  validity  of  related  circuit  calculations  to 

"small  signal")  then  the  higher  order  Taylor  series  terms,  (i.  -  I.)  , 

J    J 

become  sufficiently  small  so  as  to  become  negligible. 
Several  changes  of  variables,  such  as: 


(v  y  ■  e.i  (Ii>  V 


+  Ae 


yields  the  equations : 
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Ae. 

1 


3ei 


t  Ai-,  +  T^ 
I,  =  I,   1   oe2 


h** 


i„  =  I, 


Ae^  = 


^e2 

di 


332 

i_, 

=  L 

Ail 

+  sr 

1 

1 

2 

*2 

=  I2 

1,  = 


h«* 


i«  =  I, 


in  which  the  quantities  e.'s  and  i.'s  are  incremental  by  definition 

with  one  exception;  the  case  where  the  basic  function  is  sufficiently 

linear  for  all  the  higher  partial  derivatives  to  be  at  least  close 

to  zero. 

de. 

The  partial  derivatives  -r-r^  each  have  the  dimensions  of 

impedance.  With  this  in  mind,  the  incremental  equations  may  be  written 
in  matrix  form  where  the  definitions  of  Z.  .  are  obvious  by  inspection. 


Ae, 


Ae, 


Zll      Z12 


Z21      Z22 


"A1l" 

* 

_A12_ 

Now  we  arrive  at  the  fortunate  result  where  it  is  not  neces- 

1k 
sary  to  assume   small  signal  conditions  in  order  to  obtain  a  valid 

linear  model  of  a  transistor  using  "small  signal  Z  parameters". 

However,  it  is  necessary  to  locate  those  portions  of  the  operating 


12 


regions  for  each  transistor  where  all  four  z.  's  are  sufficiently 

linear  before  this  is  a  valid  model. 

15 
There  are  four   other  sets  of  "small  signal"  parameters 


A)  The  Y,  or  admittance  parameters 


yll     yi2 


y21     y22 


B)  The  H,  or  hybrid  parameters 


h21 


~\      — 

L12_ 

■* 

Xl 

L22 

S2 

C)  The  A,  or  transfer  parameters: 
A       B 


D)  The  M,  or  inverse  hybrid  parameters 


m 


11 


m 


m 


12 


m. 


21      22 

Each  of  the  above  coefficient  matrices,  in  its  linear 
region  may  be  arranged  to  represent  common  emitter,  common  base,  or 


13 


common  collector  orientations  through  a  series  of  well-known  trans- 

-1  c 

formations   .   In  this  respect  this  matrix  representation  is  very 
useful  for  model  construction. 

In  summary,  the  so-called  small  signal  parameter  models 
may  form  useful  quiescent  models  in  the  regions  where  all  four 
coefficients  of  each  coefficient  matrix  are  linear.  However,  lacking 
sufficient  data  on  which  to  justify  their  use,  the  systems  analyst 
should  use  a  distinctly  nonlinear  model  such  as  a  Piecewise  Linear 
Model,  Ideal  Diode  Model,  or  Regression  Model  which  is  valid  over  the 
range  of  interest. 

2.3  Piecewise  Linear  Model  Construction 

Another  class  of  models  may  be  generated  using  piecewise 

17 
linear  techniques   .   The  Piecewise  Linear  Models  are  also  valid  for  a 

wide  range  of  operating  points  provided  a  sufficient  number  of  segments 
are  employed.  An  example  of  a  Piecewise  Linear  Model  of  a  forward 
base-emitter  diode  coded  in  the  Jovial  subset  of  ALGOL  is  shown  in 
Figure  k.      In  a  computer  program,  such  a  model  consists  of  logical 
statements  defining  the  values  at  which  the  individual  segments  should 
be  switched  in  and  out  of  the  model's  network. 

Construction  of  a  Piecewise  Linear  Model  begins  with  a  graph 
of  the  objective  function  (Figure  5).  A  series  of  straight  line  seg- 
ments are  then  constructed  which,  when  appropriately  summed,  form  the 
required  objective  function  within  the  desired  accuracy.  An  example 

Piecewise  Linear  Model  for  the  I_,  versus  V^^  input  characteristic  of 

B         BE 

an  KPN  switching  transistor  is  presented  in  Figure  5B.   Several  other 
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*  1.00  FLOAT; 

*  2.00  READ  E17; 

*  3-00  Q5RS  =  0.60; 

*  4.00  ZO  =  40; 

*  5.00  FORMAT  HD5,S20,C*CHARACTERISTICS  OF  TRANSISTOR  Q5*//> 

*  6.00  FORMAT  HD5U,C*    IB     IC     HFE     RC     VSAT 

*  7 . 00      VCE     VOUT     ES */ ; 

*  8.00  PRINT  HD5; 

*  9-00  PRINT  HD5U; 

*  10.00  FOR  IB 5  =  0.005,  0.005,  0.060; 

*  11.00  BEGIN 

*  12.00  IC5MAX  =  ((E17/ZO)-0.45)/l3-0952U; 

*  13-00  IF  IB 5  GQ  0  AND  IB 5  LS  0.007; 

*  14.00  IC5  =  52.5*IB5; 

*  15.00  IF  IB 5  GQ  0.007  AND  IB 5  LS  0.0l4; 

*  16.00  IC5  =  0.18  +  26.66667*  135; 

*  17.00  IF  IB 5  GQ  0.014  AND  IB 5  LS  0.0325; 

*  18.00  IC5  =  0.3  +  17 -71429*185; 

*  19.00  IF  IB 5  GQ  0.0325  and  IB 5  LQ  IB5MAX; 

*  20.00  IC5  =  0.45  +  13. 09524* IB5; 

*  21.00  Zl  =  (E17/IC5)-Z0; 

*  22.00  HFE  =  IC5/IB5; 

*  23.00  VSAT   =  IC5*Q5RS; 

*  24.00  VCE  =  IC5*Z1; 

*  25.00  VOUT  =  E17   -  IC5*(Z1+Q5RS); 

*  26.00  ES   =  VSAT  +  VCE  +  VOUT; 

*  27.00  FORMAT  Q,F7.4,F10-3,   F9.1,2F9- s,Sf8.2,F9.2; 

*  28.00  PRINT  Q,IB5,IC5,HFE,Z1,VSAT,VCE,V0UT,ES; 

*  29.00  END 

*  30.00  FORMAT  QS,C*  */>H8,F12.4,C*  *// ; 

*  31-00  PRINT  Q£,8H(IB5MAX  =),IB5MAX; 

*  PRINT   COMPLETE 


Figure  4.   Computer  Program  for  Calculating  I  =  f(lD) 


6*  ea  70  72  74  76  78  .80 

BASE    TO    EMITTER    VOLTAGE     IV,,)-  VOLTS 

A)   C3  Emitter  Diode  Input  Characteristic 
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5       40 


SEGMENT    A  — .  _  >' 

^1-—  ^ 

** -^ i * 


0  .66  61  70  72  74  76  78  .80 

BASE    TO  EMITTER   VOLTAGE    (  V„)-  VOLTS 

B)   Piecewise  Linear  Approximation 


Segment 
A 


Equation 

v   =  0.680  +  192.85  *  r. 

BEA  B 


B 


vBE  =  0.727  +  591-39  *  iB 


v   =  0.755  +  1569.23  *  I 

■Dtp  b 


D 


^BE  =  °'7lk  +   3153*81+  *  ZB 


C)  Piecewise  Linear  Equations 


Figure  5«  Piecewise  Linear  Approximation  of  a  Function 
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1  o 

approaches  for  generating  piecewise  linear  models  are  given  in  Anner  , 

which  is  the  most  comprehensive  reference  on  the  subject.   One  con- 
straint on  the  use  of  Piecewise  Linear  Models  is  that  they  must  be 
entered  into  a  computer  program  in  segments  whose  equations  must  be 
switched  in  or  out  of  the  model  as  a  function  of  V^^  by  conditional 
program  logic  statements,  such  as  the  "IF  statement"  of  FORTRAN  Iv   , 
in  order  for  the  total  functional  value  3L  =  f(V  )  to  be  correct. 

2.h     Ideal  Diode  Model  Construction 

Suppose  a  problem  is  given  whose  solution  requires  rather 
exact  simulation  of  a  transistor's  output  characteristic  such  as  is 
depicted  in  Figure  6.   This  set  of  curves  shows  a  set  of  measured, 
static  characteristics  of  an  Emitter-Coupled-Logic  nanosecond  switch- 
ing transistor.   Contrasting  these  with  the  voltage-current  relation- 
ships of  the  controlled  sources  of  Figure  3  suggests  that  with  the 

20 
application  of  appropriate  mathematical  nonnegativity   restrictions 

one  could  form  a  suitable  model. 

An  effective  way  of  realizing  the  nonnegativity  restrictions 
is  by  adding  ideal  junction  diode  elements  to  the  model.  Thus,  by 
employing  the  diode  model  in  appropriate  series  and  parallel  config- 
urations, the  controlled  sources  for  quadrant  V-I  characteristic  may 
be  transformed  into  the  desired  form.   One  possible  algorithm  using 
the  current  and  voltage  conventions  of  Figure  7A  is : 

A)  Shunt  the  controlled  current  source  with  an  ideal 
resistive  element,  R_,  to  provide  a  modified  set  of 
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V-I  characteristics  which  have  a  linearly  increasing 
positive  slope  of  l/R  mhos  and  whose  ordinate  values 
at  V       =  0  are  unaltered  as  shown  in  Figure  7B.  Thus: 


=  1   +1 
C    1    R 


where  i  is  the  resistive  current  having  a  linear 
K 

dependency  upon  the  output  voltage  V—  since 

CB 


i„  =  V  /R_. 
R    CB  C 


B)  Inhibit  the  output  current  i_  by  placing  an  ideal 
junction  diode  element  in  series  with  the  input  branch 

such  that  current  i  now  flows  through  the  ideal  junc- 

E 

tion  diode  unidirectionally.      The  output   current  i 
becomes : 


hs  • 3  *  4 +  *B 


and  the  quasi-model  has  the  voltage  current  character- 
istic of  Figure  JC,  where  S  is  the  current  gain  of  the 
model . 

C)  Shunt  the  terminals,  across  which  the  potential 

V,_  appears  by  placing  the  ideal  junction  element  so 
CB 
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A)     B_^ 


a'c 


■*■  V 


CE 


**  (T)ai    |R 


B)      B 


E  >nC 


'CE 


c) 


D)     X 


'CE 


'CE 


•NOTE  :  VCE  =  VC,-V„ 


Figure  7.   Step-wise  Development  of  an  NPN 
Common  Emitter  Transistor  Model 
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as  to  inhibit  VnT.   from  becoming  negative.      The   output 
current  is  by  KCL : 


iC  ■  Vol  («  -^^-1)   ♦  V^R,  -  I02  (e-^/^-D 


while  the  input  current  is 


E  ""  01  K  J 


By  KCL  and  KVL  these  equations  may  be   couched  in  normal 

common  emitter  variables;   !_,  instead  of  I     and  V 

d  E      CE 

instead  of  V_.  with  the  result : 
OB 


1„  -  °F  L  +  XC0  (e  -1V°S  -   Win*.!) 


t        f-i^s,  \    t    ~qVBE/Mkt  nN        _        /    -q(    CE-  BE)Mkt  .  v 

h  =  ~hs  (1"aF)  (e  -1}  +  ^s  (e  "^ 

-vcb/rc 


which  are  the  equations  of  the  curves  in  Figure  7D. 

21 
These  relations,  according  to  Gibbons   ,  provide  a  represen- 
tation of  transistor  output  characteristics  which  ds  sufficiently  accur- 
ate that  it  falls  within  the  typical  manufacturing  spread  for  a  given 
transistor  type  over  most  of  the  normal  active  region.  Also,  one  finds 


21 


22 
upon  comparing  the  Ebers  Moll  Model  equations   with  those  of  paragraph 

C  above,  that  there  is  good  agreement  in  the  forward  active  region. 

The  above  intuitively  derived  nonlinear  model  is  easily 
extended  to  include  inverse  operation  by  adding  an  QLi^p  generator, 
as  suggested  by  the  Ebers  Moll  Model,  plus  a  shunting  ideal  resistor 
R  to  facilitate  modeling  of  the  reverse  characteristic  slope,  of 
physical  transistors,  as  shown  in  the  third  quadrant  of  Figure  6. 

The  quantities,  I  .  ,  I._,  OL  and  R  ,  in  the  above  equations 

are  variables  and,  therefore,  need  to  be  evaluated  experimentally. 

23 
Malmberg   states  that  evaluation  of  these  variables  requires  one  to 

make  a  substantial  number  of  measurements. 

2.5  The  Regression  Model 

The  Regression  Model  combines  the  simplicity  of  the  Hybrid 
parameter  model  with  the  nonlinear  modeling  flexibility  of  the  non- 
linear model  utilizing  the  ideal  junction  diode  as  discussed  in  the 
previous  section.  No  longer  is  one  restricted  to  an  incremental 
variation  about  the  "static"  or  quiescent  operating  point  as  is  true 
with  the  unmodified  Hybrid  parameter  model  of  Figure  8.   The  operating 

range  is  extended  by  equations  for  BETA  (I  )  and  V   (I  ) . 

C       BE    C 

The  first  step  in  the  modification  process  is  to  convert 

the  Hybrid  parameter  model  of  Figure  6A  to  the  ECAP  direct  current 

model  of  Figure  8.   This  is  accomplished  by  'substituting  a  battery  V^ 

B 

for  the  h      «i     voltage  generator  and  substituting  the   direct   current 
re     c  to 

gain,   BETA,   for  the   small   signal  forward  current  gain,   h„   • 
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B  o- 

+t 

Vbe 


'b 


E  o 


,hrejc 


OE 


A.  Hybrid  Parameter  Model 


B.  ECAP  Direct  Current  Model 


Figure  8.   Comparison  of  Hybrid  and  ECAP  Transistor  Models 
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The  final  step  involves  the  use  of  modeling  equations  devel- 

2k 
oped  by  Guth   "but  in  normalized  form  so  that  the  natural  logarithm 

is  taken  of  a  dimensionless  quantity.   These  modeling  equations  con- 
trol the  base  to  emitter  voltage  drop  V-^  and  th?  current  gain  BETA 
of  the  ECAP  model.   The  advantage  of  this  approach  is  that  the  cir- 
cuit analyst  need  not  use  the  ideal  junction  diode  as  a  modeling 
element.   Instead  of  the  exponential  function  and  its  attendant 

numerical  evaluation  problems  when  used  in  matrix  equations,  the  user 

25 

may  attain  the  equivalent   circuit  simulation  range  from  curves 

derived  from  either  measured  or  manufacturer's  data  points.   The 
data  points  are  fed  into  the  Linear  Regression  curve  fitting  program 
of  Appendix  C  for  the  purpose  of  generating  in  closed  form,  a  mathe- 
matical equation  which  describes  the  actual  variations  in  vL^  and 

BE 

BETA  as  a  function  of  operating  point  (collector  current,  I  )  for  each 
transistor  used.   In  this  manner  the  nonlinear  volt age -current  rela- 
tions of  any  transistor  may,  in  general,  be  established  for  any 
desired  operating  range  below  saturation. 

Note  that  the  nonlinear  Regression  Model  of  Figure  9  is  an 
extension  of  Guth's  model  for  the  purpose  of  permitting  simulation  of 
transistor  circuit  operation  in  both  the  active  and  cutoff  regions  as 
is  required,  for  example,  for  modeling  nonsaturating  Emitter  Coupled 
Logic  circuitry.   The  modification  consisted  of  adding  the  program 
logic  required  to  switch  to  a  zero  valued  function  of  collector  current, 
I  ,  whenever  the  base  current,  I  ,  decreased  below  a  preselected  refer- 
ence  level,  I  ,  or  reversed  sign  from  the  normal  forward  direction  of 
positive  amplification.  Also,  the  base  circuit  is  modified  by  setting 
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COLLECTOR 
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Ib     hie 

BASEO-t-* — WV — |  |— 

T+  +  '  '  - 

I  V 

Vbe 


hoel        (*)BETA*I 
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EMITTER  O 
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'CE 


A)  ECAP  Direct  Current  Model 


BETA  =  B  +  B.  (in  S)   +  BQ  (in  r^f 


Cx2 


vbe  =  vi +  V2  (^n  r}  +  v3  <in  r> 

0  0 


B)  Regression  Equations 


Figure  9-  Regression  Model  With  Equations 
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this  circuit  equal  to  a  branch  with  V™  =  V  ,  where  V  is  the  cutin 
^  BE    y'  7 

voltage,  in  series  with  a  resistance  of  five  megohms,  whenever 


igSo. 


Thus  the  solution  set  of  the  Regression  Model  is  now  restricted 


to  positive  values  of  BETA,  I  ,  and  I  ,  by  the  addition  of  a  conceptual 
transistor  base-emitter  diode  which  behaves  in  the  cutoff  mode  in  a 
similar  manner  to  the  ideal  junction  diode  of  the  Ideal  Diode  Model 
constructed  in  Section  2.k   above. 

The  application  of  this  Regression  Model  requires  the  user 
to  calculate  BETA  using  the  equation: 


BETA  =  T^ =£L 

B2  "  Bl 


where  BETA  is  defined  as  the  spacing   between  the  common  emitter  out- 
put characteristic  curves  of  the  Ebers  Moll  Model  operating  in  the 
normal  active  region.   There  is  an  implied  extension  of  this  defini- 
tion of  BETA  required  in  the  use  of  this  model.  For  example,  to 
extend  the  Regression  Model's  validity  into  the  cutoff  region,  it  was 

necessary  to  obtain  BETA  and  v_,„  data  at  low  currents  where  the  output 

BE 

characteristic  curves  are  highly  nonlinear.  However,  the  Regression 
Model  specifically  allows  BETA  to  be  a  variable  whereas  in  the  Ebers 
Moll  Model  BETA  is  defined  as  a  constant.   The  use  of  the  above  equa- 
tion is  justified  for  currents  which  are  near  zero,  since  the  collector 
current,  I  ,  approaches  zero  at  a  much  greater  rate  than  does  the  base 
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current,  I  .   In  the  limit  the  above  equation  does  not  converge  to  zero 

13 

as  desired,  but  it  is  a  simple  matter  to  define  BETA  as  being  zero  at 
1=0  and  1—   =   0,  since  this  is  consistent  with  all  physical  observa- 
tions of  BETA'S  behavior. 

The  base  to  emitter  forward  diode  drop,  V^,  voltage  data 

BB 

is  obtained  from  the  common  emitter  input  characteristic  curves.   Thus, 

the  user  supplies  as  input  to  the  Linear  Regression  program  a  table  of 

I_,  BETA,  and  V   which  contains  sufficient  data  to  completely  define 
C  BE 

each  transistor  in  the  region  to  be  modeled. 

The  accuracy  with  which  the  Regression  Model,  in  conjunc- 
tion with  the  ECAP  program,  will  simulate  circuit  operation  is  directly 
affected  by  the  "goodness  of  fit"  of  the  regression  equations  to  the 
experimental  data.  A  plot  of  typical  measured  physical  variations 

of  V         and  BETA  as  a  function  of  the  collector  current,  I  ,  for 
BE  0 

the  comparison  with  the  Linear  Regression  prediction  equations,  is 
shown  in  Figures  10  and  11,  respectively.   The  Linear  Regression  pro- 
gram of  Appendix  C  was  used  to  derive  these  equations.  A  convenient 
program  to  test  the  Regression  Model's  effectiveness  is  given  in 
Appendix  D. 

In  conclusion,  this  discussion  of  models  has  shown  the 
relationship  of  the  Regression  Model  to  other  types  of  transistor 
models  currently  in  wide  use.   It  is  believed  that  this  model  will 
form  the  basis  of  a  useful  parameter  variation  simulation  program. 
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3-   COMPARISON  OF  TWO 
COMPUTER  PROGRAM  CIRCUIT  ANALYSIS  METHODS 


3*1  Discussion  of  a  Basic  Circuit  Analysis  Problem 

A  problem  often  circumvented  in  the  analysis  of  transistor 
circuits  is  the  current  convergence  problem.   This  problem  arises  from 
the  appearance  of  a  dependent  generator  in  the  network  model  and  the 
simultaneous  need  to  satisfy  the  KCL  and  KVL  equations.   Specifically, 
whenever  the  circuit  external  to  the  circuit  model  has  a  direct  influ- 
ence upon  the  current  gain,  BETA,  of  one  or  more  states,  then  the 
proper  value  of  BETA  must  be  determined  prior  to  the  completion  of 
circuit  solution.   Twp  approaches  for  solving  this  problem  are:   the 
Matrix  Partition  Method  and  the  Regression  Model  Method. 

Both  of  the  above  methods  utilize  the  matrix  circuit  analysis 

27 
procedures  as  discussed  by  Seshu  .   For  ease  of  notation  interpreta- 

pO 

tion  it  is  suggested  that  one  refer  to  the  matrix  analysis   of  a 
simple  circuit  which  is  presented  in  Appendix  A. 

3 -2     Solution  of  the  Current  Convergence  Problem  -  Matrix  Partition 
Method 

The  circuit  of  Figures  12A  and  12B  is  a  proportional  con- 
trolled oven  heater  driving  amplifier  used  for  temperature  stabiliza- 
tion of  a  crystal  controlled  computer  timing  oscillator.  Using  the 
methods  of  Appendix  A,  it  is  possible  to  solve  for  all  branch  and 
element  currents. 
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The  matrix  equations  are: 


3L  =  I..  -  JG 


=  IM- 


=  AT*L„  -  JG 
M 


aT*V*m  -  JG 


=  [AT*YM*A]  *  E  -  JG 


■  tA^VA]  *  [EB  *  V  "  JG 

=  [AT*YM*A]  *  EB  +  [AT*YMM]  *  E^  -  JG 

-  [AT*YM*A]  *  Eg  +  [AT*Y *A]  *  Z*JG  -  JG 


The  dependent  collector  current  generator  matrix,  JG,  and  the  indepen- 
dent voltage  source  matrix,  E_,  are  shown  in  Figures  13  and  lk,   respec- 
tively.  Note  that  the  elements  of  the  JG  matrix  are  factored  to  show 
the  dependent  relationship  between  the  elements  of  this  matrix  and  the 

i  .  elements  of  the  I  matrix  for  which  the  solution  is  being  obtained. 
Bj  ±5 

Therein  is  the  current  convergence  problem,  referred  to  above,  in  that 
certain  of  the  unknown  currents  appear  on  both  sides  of  the  matrix 
equations.  However,  through  the  use  of  matrix  partitioning,  the  depen- 
dent current  values,  I  ,  may  be  obtained  and  the  remaining  independent 
current  values,  ]_.,  calculated.  The  procedure  is: 

Let     J  =  Q  *  In      where     Q  =  [^  k  x  1 

and    q^  =  current  gain  of  the  n   transistor,  (BETA)  . 
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Figure  13*  Dependent  Current  Generator  Matrix 
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Figure  Ik.      Independent  Voltage  Source  Matrix 
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Thus  given  the  A  and  Z  matrices  in  Figures  15  and  16  one  may  solve  for 
the  branch  current  matrix,  I  ,  where: 


IB  =  [AT*YM*A]  *EB  +  [AT-YM-A]  *Z*Q*IM  -  Q*!^ 


By  partitioning  the  I_  matrix  into  dependent   and  independent  parts 
one  obtains: 


h 


h 


d 


=     Y* 


$ 


+     Y* 


Z*Q*I 


Bd 


<S> 


Q*I 


Bd 


<S> 


where  0  is  an  appropriately  dimensioned  null  vector  in  each  case. 
Next  the  network  admittance  of  Y  matrix,  where 


Y  =  AT*Y  *A 

M 


may  be  partitioned  into  four  parts 


Y1,  Y2,  Y3,    and  Y^ 

having  the  respective  dimensions  of: 

3x3,  3  x  15,  15  x  3,  and  15  x  15 

for  the  18  branch  equivalent   circuit  of  Figure  12-B.      The  branch   cur- 
rent matrix,   upon  partitioning  the  Y  matrix,  becomes: 


h" 


** 


Bi 


Y       Y 


Y3    Yh 


$_ 
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Y        Y 

i1      x2 


Y3    \ 


Z-Q-IBc 


0 


Q*I 


Bd 


<f> 
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From  this  latter  equation  one  may  obtain  the  relations  for  I 


and  IB.. 


^d  -  VEB  +  Yl  *[Z^*V  "  ^ 


Bd 


or       IBd  =  [U  -  Y.j*Z*Q  +  Q]'1  *[Y2*Eg] 


^d    ^i  =  V*B  +  V2***^' 

where  U  is  the  identity  matrix  having  ones  at  each  position  along  the 
major  diagonal.   Two  points  should  he  noted.  First,  it  is  necessary 
to  develop  a  permutation  matrix  for  transforming  the  original  matrices 
into  the  desired  partitive  form  and,  second,  the  user  must  supply  a.' 
value  for  current  gain,  BETA,  of  each  transistor.  From  the  numerical 
analysis  viewpoint,  there  is  no  formal  algorithm  for  deriving  the  per- 
mutation matrix;  however,  through  the  use  of  cyclical  iterative  analysis 

29 

methods  '  the  exact  value  of  BETA  for  each  transistor  may  be  woven 

into  the  solution  provided  one  has  an  equation  relating  the  current 
gain,  BETA,  to  the  base  of  collector  current  of  each  transistor. 


3.3  Solution  of  the  Current  Convergence  Problem  -  Regression  Model 
Method 

In  contrast  to  the  solution  of  the  previous  section  for  the 

network  branch  currents,  the  Regression  Model  Method  utilizes  an  exist- 

30 

ing  computer  circuit  analysis  program,  ECAP  ,  which  first  solves  for 

the  network  node  voltages,  and  second,  solves  for  the  resulting  element 
voltages  and  currents. 


ko 


Mathematically  speaking  both  methods  utilize  the  same  matrix 
formulation  in  that  the  A  matrix  approach  is  used  for  converting  the 
branch  currents  into  a  statement  of  the  total  currents  entering  each 
node  as  required  by  Kirchoff's  Current  Law  (KCL). 

The  basic  cyclical  iterative  analysis  procedure  developed 

31 

by  Guth   in  conjunction  with  the  IBM  1620  version  of  ECAP  was  not 

changed  upon  coupling  the  modified  iterative  subroutine  of  Appendix 

32 
B  into  the  new  IBM  36O  version  of  ECAP  .  As  an  aid  to  understanding 

mathematical  processes  used,  the  matrix  equations  for  the  test  cir- 
cuit of  Figure  17  will  be  developed  to  show  where  Girth' s  cyclical 
iterative  technique  enters  the  circuit  analysis  computational  process. 

The  matrix  method  of  formulating  circuit  analysis  uses  a 

•  33 
fundamental  building  block  called  the  branch  element.   The  ECAP 

standard  branch  element  is  sketched  in  Figure  18.  Note  that  this 

branch  contains  provisions  for  all  possible  driving  functions  as  well 

as  the  corresponding  response  functions.   The  branch  terminal  voltage, 

V ,    is  given  by  the  matrix  equation: 
B 

VB  "  VE  "  E 


where  V  and  E  are  the  resistive  element  voltage  and  the  voltage  source 
K 

column  matrices,  respectively. 

The  branch  currents,  I  ,  flowing  out  of  the  branches,  with 
the  positive  direction  defined  as  being  from  the  m   node  to  the  n 
node  are  given  by  the  matrix  equation: 


h  -  \ +  J 


1+1 


e2=l8v 


Figure  17 .  NLAP  Common  Emitter  Test  Circuit 
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Figure  18.   ECAP  Standard  Branch 


h3 


where  I  is  the  current  through  the  branch  admittance  element,  Y. . , 
R  ij 

of  the  circuit  admittance  matrix,  Y,  and  J  is  the  sum  of  the  quantities 
J  and  J  ,  which  are  the  branch  dependent  and  independent  current  gen- 
erator matrices,  respectively. 

The  network  driving  voltages,  E,  may  be  replaced,  using 
Norton's  Theorem,  by  the  corresponding  driving  currents.  Thus  if  E 
is  the  voltage  source  column  matrix,  the  transformed  source  current 
is : 

J'  =  Y  *  E 

and  the  total  driving  current  in  the  network  becomes: 


JT  =  JI  +  JD  +  Y*E* 


Kirchoff's  Current  Law  in  matrix  form  is 


I  =  A*I  =  0 
N      B 


where  now         I£  =  IR  +  Jy, 


and  thus,  by  expansion: 


IN  =  A  *  IB  =  [A  *  IR]  +  [A  *  JT] 


Consequently:      A  *  Y  *  V  =  [-A]  *  J_ . 

K  1 


From  Kirchoff's  Voltage  Law  in  matrix  form: 


VE  "  ^  *  V 


hk 


Combining  the  two  previous  equations 


[A  *  Y  *  A  ]  *  VN  =  [-A]  *  J  j 


and  solving  for  the  node  voltage  matrix,  V, 


N 


VN  =  [A  *  Y  *  AT]_1  *  [-A]  *  [JT] 


=  [A  *  Y  *  A  J"1  *  [-A]  *  [J  +  J  +  Y  *  E] 


Following  through  now  with  the  illustrative  example,  we  begin 
by  writing  the  KCL  equations  for  the  circuit  of  Figure  17. 


0 


-i2     -i3     -iu      0        0 


0 


\  ° 


0        0         0        +1.  m  0 


0      0 


-i,     +±c         -is       =   0 
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These  equations  become  in  matrix  notation  I  =  A  *  I  ,  or 
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_Inl" 

zm 
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and: 


h5 


The  voltage  and  current  source  matrices  are,  respectively: 


- 

el 

e2 

e3 

% 

e5 

e6 

_     - 

18.0 

18.0 
o 


V. 


JT   =  JZ  +  JD  +  ^E 


18.0 
18.0 

0 
-1.25 

0 

0 


$     + 


Jm  =  $   + 


JD1 
^D2 
JD3 

^D5 
JD6 


0 
0 
0 
0 
0 
BETA  *  i 


Bk 


[y*e] 


r 


+  [Yj 


18.0 
18.0 

o 
-1.25 

0 

0 


1*6 


Note  that  the  voltage  source  matrix,  E,  clearly  shows  the 

role  of  V--,  and  that  the  dependent  current  generator  matrix,  J, 
BE  D 

clearly  shows  the  role  of  the  transistor's  current  gain,  BETA.   These 
are  the  two  quantities  which  are  iteratively  changed  by  the  Regression 
Model-  However,  it  is  very  difficult  to  solve  the  matrix  equation, 


VN  =  [A*Y*AT]  _1  *   [-A]   *  [Jz   +  JD  +  Y*E], 


when  the  dependent  current  BETA*:L>  is  not  known.  For  this  reason, 
it  is  necessary  to  modify  the  above  equation  slightly. 

The  requisite  modification  is  obtained  by  writing  the  node 
equations  for  the  circuit  of  Figure  17 •  Upon  substituting  the  equiva- 
lent expression  for  i  ,  in  terms  of  V  ,  V  ,  V  ,  and  R,  and  collecting 
terms  one  obtains  the  matrix  relation: 


Y  *  V  =  J 
N    N    T 


3^ 
where  Y„  is  the  nodel  conductance  matrix   of  Figure  19  and  satisfies 

the  equation: 


Y  =  [A*Y*AT] 


which  was  discussed  above.   The  corresponding  nodal  equivalent  current 
vector  is  given  in  Figure  20.   This  current ' vector  satisfies  the 
relation: 


JT  "  JI  +  JD  +  Y*E 


which  was  derived  above. 


G 
O 


hi 


G 

O 
O 

ir\ 

H 


G 

~o 

H 

+ 

G 

O 
O 

I 


H 

+ 

G 

O 

C— 
OJ 

+ 

G 

O 
O 


H 

+ 

G 

O 
O 

m 


o 

CQ 

>> 

H 

-P 

•H 
O 

•H 

H 
ft 
X 


o 
o 


G 

O 
O 
O 


G 

H 

H 
+ 
G 
O 


o 
o 


G 


^-^ 

o 

+ 

G 

o 

LT\ 

G 

O 

^-^ 

O 

H 

O 

Lf\ 

O 

■* — ^ 

+ 

LTN 

< 

H 

G 

~H~ 

H 

pq 

O 

+ 

O 

G 

<~ 

o 

EH 

O 

W 

r- 

pq 

j- 

i 

~H~ 

1 

i 

pq 

^S 
-P 
•H 

■P 


o 

-p 

o 

I 

o 
o 

■d 

o 


ON 
H 

•H 


>H 


kQ 


! 

a    ' 

o 

o 

IT\ 

<~ 

G 

G 

1 

O 

O 

pq 

O 

o 

ir\ 

LTN 

* 

W 

^H 

w 

pq 

pq 

pq 

> 

> 

> 

G 

O 
O 

o 

CM 
H 

00" 
H 


G 
O 

m 
en 

GO 


o 
o 


pq 


> 
i 


h> 


k9 


3-U  Cutoff  Conditions  and  the  Regression  Model 

The  conditions  which  must  be  satisfied  for  each  transistor 

in  a  network,  which  is  in  a  nonconducting  or  cutoff  state,  are  BETA 

equal  to  zero  and  a  base  to  emitter  voltage,  V ,  less  than  or  equal 

BE 

35 
to  the  cutin  voltage   ,  V  •   The  condition  of  BETA  =  0  implies  that  the 

collector  current,  I  ,  is  zero,  provided  the  collector  leakage  amount 

I   ,  is  negligible,  since  at  low  values  of  I  : 

Total  I 

BETA  =  Total  I 
B 

This  further  means  that,  when  there  is  no  significant  base  current 

flowing,  V   is  determined  by  the  circuits  external  to  each  transistor. 
BE 

From  Figure  9,   the  Regression  Model  equations  are: 


BETA  =  Bn  +  B0  (in  ^  )  +  B   (in  -£  )2 
1  X0      3     ^ 


VBE   =  Vx  +  V2  (in  j£  )  +  V3  (in  ^  f 


Since  the  ECAP  model  is  a  strictly  linear  network,  the  currents  in 
the  base,  emitter  and  collector  circuit  branches  may  flow  in  either 
the  conventional  positive  direction  for  an  NPN  transistor  or  in  the 
opposite  direction. 

A  transistor's  current  gain,  BETA,  may  become  zero  in  one  of 
three  ways.  First,  BETA  may  be  made  zero  by  using  a  logical  operator 
statement,  such  as  the  FORTRAN  IV  code: 


50 


IF  (I    .LE.   CUTOFF)  BETA  =  0.0 

placed  in  Subroutine  PIRATE  after  the  BETA  Regression  Model  equation 
has  been  evaluated.   Second,  it  may  become  zero  by  carefully  selecting 
the  BETA  regression  equation  which  predicts  a  zero  value  for  BETA  at 
a  suitably  low  value  for  the  collector  current,  I  .   Third,  a  combina- 
tion of  the  previous  two  methods  may  be  used.   The  first  method  tends 
to  generate  discontinuities  which  make  convergence  difficult  whenever 
CUTOFF  is  an  arbitrarily  selected  current  level.   The  second  method 
gives  a  smooth  curve  down  to  very  large  negative  values  of  BETA,  as 
I  becomes  increasingly  smaller.  An  example  of  this  event  is  shown  in 
Table  1.   Note  that  in  the  BETA  equation  of  Table  1,  CUTOFF  corresponds 
to  a  value  of  I  somewhere  between  1.0  and  12.0  microamperes  and  I_ 
corresponds  to  1.0  microampere.   The  third  method  involves  a  combination 
of  the  simplicity  of  the  previous  two,  with  CUTOFF  chosen  from  the  BETA 
equation  and  with  the  "IF  statement"  inhibiting  negative  values  of  BETA. 

A  transistor's  base  to  emitter  voltage,  VRT?,  is  a  function 
of  I  throughout  the  active  region.   The  quality  of  fit  for  a  typical 
Regression  Model  V-^  equation  is  shown  in  Table  2.  Note  that  with  the 

same  value  of  CUTOFF  used  for  calculating  V  ,  one  obtains  a  number 

BE 

close  to  V  .  However,  as  the  drive  to  any  transistor  falls  below 

V__  =  V  ,  the  voltage  across  its  emitter  to  base  terminals  becomes  a 
BE    7 

function  of  the  external  circuitry  surrounding  it.   Thus  at  CUTOFF 

and  below,  a  logical  switch  needs  to  be  incorporated  into  Subroutine 

ITRATE;  so  that  the  V__  may  be  a  function  of  some  external  circuit 

Bill 

current,    I    .     A  better   external   circuit  dependent  variable  would  be 
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a  voltage,  but  this  leads  to  the  matrix  solution  problem  similar  to 
that  discussed  in  Section  3-2,  where  now  node  voltage  is  the  indepen- 
dent variable.  Note  that  both  Table  1  and  Table  2  were  calculated 
using  the  Regression  Model  Evaluation  Program  listed  in  Appendix  D. 

3*5  Extension  of  the  Regression  Model 

Two  approaches  requiring  only  minor  Regression  Model  modi- 
fication were  explored.  First,  V^^   was  set  equal  to  the  cutin  voltage 
with  the  result  that  satisfactory  predictions  of  network  behavior  can 
be  made  for  operating  conditions  close  to  cutoff.   Obvious  inaccur- 
acies are  introduced  as  each  transistor's  operating  point  passed  well 
into  the  cutoff  state.   The  second  approach  consisted  of  adding  a 
large  impedance,  ZOFF,  in  series  with  the  base  circuit  of  the  transis- 
tor for  the  purpose  of  simulating  ideal  diode  behavior.  The  ZOFF 
impedance  of  5> 555* 555*0  was  switched  into  the  circuit  by  a  logical 
operation  statement  in  Subroutine  ITRATE  whenever  the  I  current 
became  less  than  I  •   This  latter  approach  was  tested  using  the  cir- 
cuit of  Figure  1-B  with  the  result  that  current  convergence  was  not 
obtainable. 

A  check  was  made  of  the  effects  of  setting  BETA  equal  to 
zero  and  V-,^  equal  to  V  ,  by  examining  the  matrices  used  in  the  matrix 
equation: 


Y  *  V  =  J  . 
N    N    T 


5h 


The  solution  to  this  set  of  simultaneous  equations  will  exist  provided 
the  Y  matrix  has  an  inverse.  By  inspection  of  Figure  19,  it  is 
apparent  that  this  will  always  be  the  case  provided  the  admittances 
are  real  and  positive;  since  BETA  appears  only  in  nonmajor  diagonal 
terms.   The  possible  range  of  variations  in  V--,  also  will  not  affect 

-Dili 

the  existance  of  Ym  inverse  since  V   only  appears  in  the  J  matrix 

i\  BE  T 

of  Figure  20.   This  means  a  solution  for  V  will  exist  even  if  J  is 
a  zero  vector. 


3«6   Conclusion 

Summarizing  briefly,  the  computer  program  MAP  solves 
for  the  linear  current  in  each  of  the  branches.   The  resulting 
collector  current  is  then  used  as  the  dependent  variable  in  the 
Regression  Model  equations  of  Figure  9  hy  Subroutine  ITPATE.   The 
iterative  process  is  started  in  accordance  with  Guth's  flow  chart 
of  Figure  21  until  convergence.   Then  the  proper  values  of  BETA  and 
V™  are  obtained  as  a  function  of  the  collector  current  in  the  active 

-D-tLi 

region.   The  calculations  made  by  the  program  may  be  checked  both 
before  iteration  starts  and  after  convergence  is  obtained  by  using 
the  PRINT  MI  input  code  described  in  4.3 •  A  set  of  matrices  obtained 
in  this  are  shown  in  Figures  22  and  23,  respectively. 

One  interesting  example  of  the  accuracy  of  NLAP  is  a  DC 
option  analysis  of  the  linear  test  circuit  of  Figure  17 .   The  NLAP 
computed  results  are  compared  with  manually  calculated  results,  for 
the  nodal  conductance  matrix,  in  Figure  2k,    and  for  the  equivalent 
current  vector  matrix  in  Figure  25 •   The  results  compare  closely  out 
to  the  sixth  decimal  place. 
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Obtain  a  nominal 
solution 


Calculate  new  values 
for  BETA  and  VBE  using  values  of 
collector  currents  just  calculated 


Yes 


Yes 


I 


Insert  new  values  of  BETA  and 
VBE  into  the  proper  matrices 


Calculate  a  new  set  of  collector 
current  values 


Calculate  the  change  in  value 
of  preselected  reference 
current 


Calculate  the  changes  between 
previous  solution  and  the 
present  solution 


END 


Figure  21.   Flow  Chart  of  Iterative  Regression 
Model  Solution  Method. 
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NODAL  CONDUCTANCE  MATRIX 
ROW  COLS 

11-4  o. 238679200-01-0. U9999990D-02  0.0        0.0 
15-5  0.0 

2  1-1+  -O.50U99988D  00  0.10050732D  03-0.13333332D-03-0.21758056D-02 

2  5-5  -0.10000002D  03 

3  1   _  4       O.I+9999988D  00-0.50013321D  00  0. 1000001 5D  03  0.0 
3     5-5       0.0  ' 

k    1  -  k     0.0       -0.21758056D-02  0.0        0.10000219D  03 
h    5  -  5   0.0 

5     1   -   k       0.0  -0.10000002D  03  0.0  0.0 

5     5-5       0.10000003D  03 


EQUIVALENT  CURRENT  VECTOR 

NODE  NO.  CURRENT 

1  0.22717916D-01 

2  -0.38884987D  00 

3  0.11538497D  03 
h  -0.29000000D  03 
5  0.0 


NODAL  IMPEDANCE  MATRIX 
ROW  COLS 


1 

1 

-  k 

1 

5 

-  5 

2 

1 

-  k 

2 

5 

-  5 

3 

1 

-  1+ 

3 

5 

-  5 

k 

1 

-  h 

h 

5 

-  5 

5 

1 

-  h 

5 

5 

-  5 

0.5293588UD  02  0.521722U8D  00  0.69562887D-06  0.11351i+l8D-04 

0.521722U3D  00 

0.52693899D  02  0.21+901+866D  01   0.332061+35D-05   0.5^1869600-04 
0.21+904861+D  01 
-0.11396619D-02  0.984712U3D-02  0.99999983D-02  0.21^2l+959D-06 
0.98471233D-02 

0.11464917D-02  0.5U186960D-04   0.72249166D-10  0.99997821D-02 
0.5U186955D-04 

0.52693894D  02  0.24904864D  01  0.33206432D-05   0.5^1869550-04 
0.25004861D  01 


Figure  22.      DC  Solution  Matrix  Obtained  Prior  to  Iteration 
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NODAL  CONDUCTANCE  MATRIX 
ROW  COLS 

li-U  0.23867920D-01-0.U9999990D-02  0.0         0.0 

15-5  0.0 

2     1-1+  -0.49911976D  00  0.10050l¥+D  03-0.13333332D-03-0.21758056D-02 

2  5-5  -0.10000002D  03 

3  1   -  h  O.49I+II976D  00-0.49425310D  00  0. 1000001 5D  03  0.0 

3  5-5  0.0 

h    1  -  k  0.0        -0.21758056D-02  0.0         0.10000219D  03 

4  5-5  0.0 

5  1  -  k  0.0  -0.10000002D  03  0.0  0.0 
5  5-5  0.10000003D  03 


EQUIVALENT  CURRENT  VECTOR 

NODE  NO-  CURRENT 

1  0.22654877D-01 

2  -0.37802932D  00 

3  0. 11537^210  03 
k  -0.29000000D  03 
5  0.0 


NODAL  IMPEDANCE  MATRIX 
ROW  COLS 


1 

1 

-  h 

1 

5 

-   5 

2 

1 

-  h 

2 

5 

-  5 

3 

1 

-  k 

3 

5 

-  5 

k 

1 

-  k 

h 

5 

-  5 

5 

1 

-  k 

5 

5 

-  5 

0.52935133D  02  0.52783296D  00  0.70377616D-06  0.11U84368D-04 

0.52783291D  00 

0.52690314D  02  0.25196555D  01  0.33595354D-05  0.5^82l605D-OU 

0.25196553D  01 

-0.1139M+18D-02  0.98453338D-02  0.99999983D-02  0.21U21063D-06 
0-98453328D-02 

0.11U64137D-02  0.5i+82l605D-Ol+  0.73095358D-10  0-99997821D-02 
0.5^82l600D-0U 
0.52690309D  02  0.25196553D  01  0.33595350D-05   0. 5i+82l600D-04 

0.25296550D  01 


Figure  23 •   NL  Solution  Matrix  Obtained  After  Completion  of  Iteration 
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N 


+0.00229606       +0.00000000      -0.00200000 


+0.10000000     +0.003030^0     -0.10000000 


■0.10200000 


+0.00000000     +0.10570380 


a)  Manual  Calculations 


"N 


+0.22960991D-02    0.00000000D+00  -0.19999999D-02 


+0.9999996+D-01     0.30304029D-02   -0.10000000D+00 


-O.IOI99996D+OO    -0.99999966D-07   +0.10570377D+00 


b)  MAP  Calculations 


Figure  2k.      Nodal  Conductance  Matrix  Calculations  for  BETA  =  50 
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+0.00400000 


+0.1795i+51+5 


-0.12750000 


a)     Manual  Calculations 


+O.39999976D-O2 


+0.1795i+527D+00 


-0.127i+9982D+00 


b)     NLAP  Calculations 


Figure  25.     Nodal  Current  Vector  Calculations 

for  BETA  =  50  and  V^  =  -1.25  Volts 

is  ill 


60 


A  second  example  of  the  accuracy  of  NLAP  is  a  E  option 
analysis,  the  output  of  which  is  compared  with  laboratory  measured 
performance  in  the  graph  of  Figure  26  for  the  emitter  follower  of 
Figure  27. 
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N, 


[OR  =  459  XI   1 
[NOR  =  1068X1 J 
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r8  y 
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•I 


Figure  27.   OR  Emitter  Follower  Schematic  Diagram 
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k.     APPLICATION  OF  NLAP  TO  LOGIC  CIRCUIT  ANALYSIS 

4.1  Introduction 

NLAP  is  intended  to  be  used  for  analyzing  nonsaturating 
logic  circuitry  such  as  emitter  coupled  logic  blocks.   The  capability 
now  exists  within  the  program  for  linear  region  parameter  variation 
analysis.   Should  the  user  desire  to  perform  special  calculations,  a 
subroutine  may  be  coded  and  called  just  prior  to  the  return  of  control 
to  the  main  program  by  Subroutine  ITPATE.  The  NLAP  calling  procedure 
and  its  associated  system  procedure  is  given  in  Figure  28. 

k.2     Preparation  of  NLAP  Input  Data  Deck 

The  input  language  to  NLAP  is  similar  to  ECAP's^  and 

is  aligned  closely  with  the  language  of  the  electronic  engineer.   Each 

of  the  circuit  elements  use  a  simple  descriptive  letter  coding  and  the 

topology  of  the  circuit  interconnections  is  identified  by  a  sequence  of 

branch  and  node  numbers.  The  polarity  of  positive  current  flow  is 

inserted  by  ordering  the  nodal  connection  using  identifying  numbers 

in  a  from  node  m  to  node  n  arrangement. 

The  general  branch,  or  B.,  statement  consists  of  a  branch 

J 

number,  j;  the  from  and  to  node  numbers,  m  and  n,  respectively;  the 
braich  resistance,  R,  in  ohms;  the  independent  voltage  source  E,  in 
volts;  and  the  independent  current  source  I,  in  amperes.   In  other 
words,  the  general  branch  statement  is: 


B.  =  N(m,n),  R=7,  E  =  e,  I  =  i 

J 
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/*FORMAT  PR,DDNAME=FT10F001,COPIES=l 


EXEC 


//JOBLIB 
//  EXEC 
//MAP 

II 
II 

'I, 

//FT05F001  DD 
//FT06F001  DD 
//FT07F001  DD 
//FT08F001  DD 
//FT09F001  DD 
//FTlOFOOl  DD 
//FTllFOOl  DD 
//FT12F001  DD 
//MAP.SYSIN 

/* 


DD     DSNAME=SYS1.CKTLIB,DISP=( OLD, PASS) 

MAP 

PGM=MAP 
* 

EXECUTE  MAP       --       JOBLIB  CARD  REQUIRED 
//JOBLIB  DD  DSNAME=SYS1.CKTLIB,DISP=( OLD, PASS) 


DD 


DDNAME=SYSIN 

SYSOUT=A,DCB=(RECFM=FBA,LRECL=133,BLKSIZE=798) 
UMET=(CTC, , DEFER) ,DCB=(RECFM=F,LRECL=80,BLKSIZE=80) 

UNIT=(CTC, , defer) ,dcb=(recfm=f,lrecl,8o=blksize=8o) 

UNIT=(CTC,, DEFER) =DCB=(RECFM=F,LRECL=80,BLKSIZE=80) 

UM!T=(CTC, , DEFER) ,dcb=(recfm=f,lrecl=8o,blksize=8o) 

UNIT=(CTC, , DEFER) ,DCB=(RECFM=F,LRECL=80,BLKSIZE=80) 

UMET=(CTC, , DEFER) ,dcb=(recfm=f,lrecl=8o,blksize=8o) 

*  MAP  DATA  DECK  GOES  fflIRE 


Figure  28.     MAP  Calling  Sequence 
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where  j,  m,  n,  y ,  e  and  i  are  place  holders  for  specific  numeric  values 

which  vary  from  "branch  to  branch.  An  example  of  the  specific  coding 

for  the  circuit  of  Figure  27  is  listed  in  Figure  29. 

The  dependent  current  sources  are  denoted  by  a  T  statement 

containing  the  controlling  and  the  controlled  branch  identification 

numbers,  wotether  with  the  associated  current  gain  value,  BETA. 

The  regression  equation  coefficients  for  the  BETA  and  V^- 

•BE 

equations  are  entered  in  T  and  B.  statements,  respectively,  which 

k      0 

are  of  the  form: 

TR     BETA  (i)  =  B1,  B2,  B3 

B.  E(i)     =  V1,  V2,  V3 

where  Tv  and  B  .  start  in  card  column  one  with  the  BETA  (p)  and  the 
k      J 

E(p)  dependent  variable  statements  starting  column  seven. 

The  cyclical  iteration  process  is  controlled  by  two  succes- 
sive statements 

CURRENT  -  7 


and 


TOLERANCE 


which  imply  that  the  iteration  process  is  to  terminate  whenever  the 
newest  solution  for  the  current  in  branch  7  changes  by  less  than  the 
prespecified  amount,  5,  from  the  previous  iterations  calculated  value. 
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ID21        OR       EMITTER  FOLLOWER,    RE  =  1+59-6  OHMS        -DCAP 

DC 
Bl         N(0,l),R=53  ,E=1.0 

B2  N(l,2),R=200.0,E=-0.770 

B3         N(0,3),R=0.01,E=+1.15 

B*+  N(U,0),R=0.01,E=+2.90 

B5         N(3,2),R=7500 

B6         N(5,^),R=^59-6 

B7         N(2,5),R=0.01 

B8         N(5.0,R=100000 

Tl  B(2,5),BETA=100 

PRINT  NV,BV,BA,BP,MI 

EX 

END 


A.      DC  Circuit  Analysis  Data  Deck 


ID21 

OR       EMITTER  FOLLOWER,   RE  =  459-6  OHMS        -NLAP  DATA 

NL 

Bl 

N(0,l),R=53                     ,E=+1.0 

B2 

N(l,2),R=200.0,E=-0.770 

B3 

N(0,3),R=0.01,E=+1.15 

Bk 

N(U,0) ,R=0.01,E=+2.90 

B5 

N(3,2),R=7500 

B6 

N(2,l+),R=l+59.6 

B7 

N(2,5),R=0.01 

B8 

N(5,0),R=100000 

Tl 

B(2,5),BETA=100 

PRINT     NV,MI 

EX 

Tl 

BETA(5)  =  -551+- 67056,  +172.1+2553,  -11. 1718 

B2 

E(5)=-0. 5^023,  -0.021+29,  -0.00021 

CURRENT  =  6 

TOLERANCE  =0.01 

EX 

RETURN 

END 

B.      NL  Circuit  Analysis  Data  Deck 


Figure  29.     Typical  NLAP  Input  Data 


67 


The  amount  of  tolerable  deviation  from  current  convergence  is  stated 
in  a  decimal,  not  percent,  format.   One  example  is: 

TOLERANCE  =0.01 

which  means  current  convergence  to  within  one  percent. 

^•3  Solution  Control  Codes 

Two  types  of  analyses  are  possible  using  NLAP,  linear  direct 
current  and  iterative  nonlinear  direct  current.   The  solution  control 
cards  are  DC  and  NL,  respectively. 

A  DC  problem  deck  must  contain  only  those  cards  up  to  the 
first  EX  or  EXECUTE  card  as  shown  in  Figure  29-A.   The  full  set  of 
ECAP/36O  DC  solution  options  '  SE,  WO,  ST,  and  MI  are  available  to  ' 
the  DC  user.   These  options   obtain  for  the  user:  element  sensitiv- 
ities, worst  case  tolerance  analysis,  standard  deviations  of  the  node 
voltages  assuming  component  values  which  are  statistically  independent, 
random,  and  normally  distributed  about  their  mean  values,  and  a  set 
of  solution  matrices,  respectively. 

An  NL  problem  deck  must  contain  the  same  network  description 

cards  as  the  equivalent  DC  problem  deck,  plus  the  regression  equation 

cards  for  BETA  and  V_„  and  the  CURRENT,  TOLERANCE,  EXECUTE,  RETURN  and 

BE 

END  cards  as  shown  in  Figure  29 -B.   These  last  three  cards  are  self 
descriptive  and  the  user  should  observe  that  the  'END  card  is  used  to 
s  ignal  the  end  of  the  input  data  so  that  control  may  be  passed  to  the 
system  monitor. 

It  should  be  noted  that  the  only  case  where  the  PR  or  print 
card  is  essential  is  when  the  user  desires  to  print  out  the  nodal 
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conductance  matrix,  the  equivalent  current  vector,  and  the  nodal 
impedance  matrix.   To  obtain  these  matrices  the  required  code  is: 

PRINT  NV,  MI 

This  coding  was  used  to  obtain  the  matrices  shown  in  Figures  22  and 
23- 

h.k     Solution  Output  and  Identification 

When  the  user  desires  to  perform  a  parameter  variation 
analysis,  careful  attention  should  be  given  to  the  preparation  of  the 
title  card  which  is  placed  at  the  beginning  of  each  problem  deck  for 
solution  identification  purposes. 

The  letters  ID  are  placed  in  card  columns  one  and  two  for 
the  purpose  of  signal  by  the  subroutine  ITRATE  to  search  for  subse- 
quent code.   Card  column  three  is  reserved  for  a  problem  identification 

number  from  one  to  nine  which  is  read  and  written  by  the  program  in  Al 

39 
format   .   Card  column  5  is  reserved  for  a  flag  which  resets  the 

solution  number  to  one  after  a  series  of  related  solutions  have  been 

run  for  a  given  problem  and  the  user  desires  to  switch  to  a  new 

solution  series  of  the  same  or  different  problem.   Any  integral 

value  between  one  and  nine  will  reset  the  solution  counter  to  one. 

There  are  a  variety  of  output  options  which  may  be  selected 

through  the  use  of  IBM  36O  Job  Control  Language   .   The  two  principal 

modes  of  output  communication  are  80  column  card  images  and  133 

column  printed  text.   In  planning  the  FTXXF001  data  sets  listed  in 

the  system  procedure  of  Figure  2k,   the  user  was  given  a  very  flexible 
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list  of  data  sets,  numbers  FT07  through  FT12,  upon  which  to  call  and 
use  for  parameter  variation  studies.  As  a  result,  these  later  data 
sets  are  all  80  bytes  long  and  are  blocked  at  80  bytes  so  that  the 
user  may  employ  either  cards  or  directed  access  devices  for  data 
transmission  to  an  appropriate  plotter.   The  MAP  output  routine 
ECB25  is  listed  in  Appendix  E.   Thus,  the  user  may  select  the  required 
data  set  by  consulting  this  subroutine  and  the  definitions  of  the 
array  names  given  in  Table  3* 

A  typical  MAP  solution  obtained  without  use  of  the  PRINT 
statement  of  ECAP  is  contained  in  Figure  30  and  is  the  output  of  the 
FT06F001  data  set. 

4.5  Conclusion 

The  range  of  validity  of  the  Regression  Model  has  been 
found  to  be  limited  to  the  active  region.  The  modification  of  these 
equations  and  the  subsequent  testing  were  outside  the  scope  of  this 
investigation  and  are,  therefore,  a  subject  for  further  study. 
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OUTPUT  VARIABLE 


NLAP  STORAGE  ARRAY 


Branch  Voltage 

Branch  Current 

Element  Voltage 

Element  Current 

Branch  Power  Dissipation 

Node  Voltage 

Transconductance 

BETA 

Element  Admittance 

Element  Impedance 


CCSAV  (1,  J) 

CCSAV  (2,  J) 

CCSAV  (3,  J) 

CCSAV  (k,  J) 

CCSAV  (5,  J) 

CCSAV  (6,  J) 

CDSAV  (1,  J) 

CDSAV  {2,  J) 

Y(J) 

X(J) 


Table  3*      Summary  of  NLAP  Output  Variable  Storage  Arrays 
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APPENDIX  A 
A  NETWORK  ANALYSIS  PROCEDURE 

1.  Introduction 

The  formulation  of  an  electronic  circuit's  mesh  solution 
begins  by  writing  the  branch  level  matrices  of  Figure  A2  for  the  net- 
work shown  in  Figure  Al.  Each  of  these  matrices  is  a  column  matrix 
containing  one  row  element  for  each  network  branch. 

2.  Matrix  Equation  Formulation 

The  impedance  matrix,  Z,  of  Figure  A3  is  generated  using 
the  relation: 


Matrix  element  Z..  =  Z  ,    ,  . 
n     branch. 

l 


The  Kirchoff's  Voltage  Law  coefficient  matrix,  A,  of  Figure 
A3  is  formed  by  summing  the  voltages,  V.,  around  each  mesh  with  the 
algebraic  sign  determined  by  the  agreement  (+),  or  disagreement  (-), 
with  the  arbitrarily  preselected  positive  current  direction.   This  A 
matrix  has  as  many  rows  as  the  network  has  meshes  and  as  many  columns 
as  the  network  has  branches. 

The  quantities,  Z,  JG,  E,  and  A,  are  processed  by  a  series 
of  matrix  equations  to  arrive  at  the  solution  for  the  network's 
branch  current  matrix,  I  .  These  equations  are: 

a)  All  Norton  equivalent  current  sources  are  converted 
Thevenin  equivalent  voltage  sources: 
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Figure  A2.  Branch  Level  E,  JG,  V,  and  I  Matrices 
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N  =  Z  *  JG 

b)  The  converted  current  sources  are  summed  with  the 
native  Thevenin  equivalent  voltage  sources : 

EE  =  E  +  N. 

c)  The  voltage  sources  around  each  mesh  are  summed  by  a 
matrix  equation: 

EM  =  A  *  EE. 

d)  The  mesh  impedance  matrix  is  formed  by  transposing  the 

A  matrix,  premultiplying  the  Z  matrix  by  A,  and  then  post- 

T 
multiplying  A  *  Z,   the  A  *  Z  product  by  A  : 


AM  =  A  *  Z  *  AT. 


e)  The  network's  mesh  current  is  calculated  by  premulti- 
plying the  EM  matrix  by  the  inverse  of  the  ZM  matrix: 

IM  =  YM  *  EM. 

f )  The  final  step  is  to  obtain  the  individual  currents 
in  each  branch  by  summing  up  the  mesh  currents  passing 
through  each  branch  and  then  subtracting  the  current 
generator  matrix,  JG: 


IB  =  AT  *  m   -  JG. 
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3«   Evaluation  of  Matrix  Equations 

The  above  matrix  equations  are  herein  evaluated  for  the 
network  of  Figure  Al. 


a) 


N  =  Z  *  JG  = 
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d) 


AM  =   [A*Z*A  ] 
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e) 

IM  = 


IM 


3 


=  YM  *  EM  = 


10 
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f) 


IB  =  A     *   IM  -  JG  = 


-10     0 

-1  +1     0 

-1.88' 

0-10 

0     0+1 

-0.397 

0  +1   -1 

-1     0  +1 

+0.036 

-  [J]  = 


+1.88 

+1.1+8 

+O.1+0 

+O.036 

-0.1+1+ 

+1.1+2 


-0.12 
+1.1+8 
+1+.60 
-1.96 
-0.1+1+ 
+1.92 


A  network  analysis   computer  program  utilizing  the  above 
matrix  equations   is  listed  in  Figure  Al+,   with  a  sample  output   shown 
in  Figure  A5« 


amps 
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*  1.00  FLOAT; 

*  2.00  ARRAY  A  3  6  F;  ARRAY  E  6  1  F; 

*  3.00  ARRAY  JG  6  1  F;  ARRAY  Z  6  6  F; 

*  4.00  ARRAY  N  6  1  F;  ARRAY  EE  6  1  F; 

*  5-00  ARRAY  EM  3  1  F;  ARRAY  AT  6  3  F; 

*  6.00  ARRAY  ZAT  6  3  F;  ARRAY  AM  3  3  F; 

*  7-00  ARRAY  YM  3  3  F;  ARRAY  IM  3  1  F; 

*  8.00  ARRAY  IBM  6  1  F;  ARRAY  IB  6  1  F; 

*  9.00  ARRAY  IBX  6  6  F;  ARRAY  IBS  6  1  F; 

*  10.00  ARRAY  PD  6  1  F;  ARRAY  P  6  6  F; 

*  11.00  ARRAY  RP  6  6  F;  ARRAY  S  6  1  F; 

*  12.00  ARRAY  SP  6  1  F; 

*  13-00  A[0,0]=-1;  A[o,iJ=-i;  A[o,5l=-i; 

*  14.00  A[l,l]=  1;  A[l,2]=-1;  A[l,4]=  1; 

*  15-00  A[2,3]=  1;   A[2,4]=-l;  A[2,5]=  1; 

*  16.00  E[l,0]=10; 

*  17.00  JG[0,0]=  2;  JG[2,0]=  5;  JG[3,0]=  2; 

*  18.00  RP[0,0]=  1;  RP[1,1]=  1;  RP[2,2]=  1; 

*  19.00  Z[0,0]=  5;  Z[l,l]=  2;  Z[2,2]=  2;  Z[3,3]=  5; 

*  20.00  A[3,3]=  5;  Z[4,4]=  5;  Z[5,5]=  h; 

*  21.00  N  =  Z  *  JG;  EE  =  E  +  N;  EM  =  A  *  EE; 

*  22.00  FOR  I  =  0,  1,  2;   BEGIN  FOR  K  =  0,  1,  5;   BEGIN 

*  23.00  AT[K,  I]  =  A[I,  K];   END   END 

*  24.00  ZAT  =  Z  *  AT;  AM  =  A  *  ZAT;  YM  =  ZM**  -1.0; 

*  25.00  IM  =  YM  *  EM;  IBM  =  AT  *  IM;  IB  =  IBM  -  JG; 

*  26.00  IBX[0,0]  =  IB[0,0]; 

*  27.00  IBX[3,3]  =  IB[3,0]; 

*  28.00  IBS  =  IBX  *  IB; 

*  29.OO  PD  =  Z  *  IBS;  P  =  RP**  -1.0; 

*  30.00  S  =  P  *  PD;  S  =  100  *S; 

*  31-00  FORMAT  LINE,  S27,  C*MATCK  DATA  OUTPUT*// 

*  32.00  PRINT  LINE; 

*  33-00  FORMAT  HDR,  C*     I     Z     IB     PD     RP 

*  34.00  S         SP*/; 

*  35-00  PRINT  HDR; 

*  36.00  FORMAT  UNITS,  C*  BRANCH  OHMS   AMPS   WATTS 

*  37-00      WATTS       NONE      PERCENT*/; 

*  38.00  PRINT  UNITS; 

*  39-00  FORMAT  MATCK,  16  ,6(Flo.2); 

*  40.00  FOR  I  =  0,1*5;   BEGIN 

*  41.00  PRINT  MATCK,   I,Z[l,l] ,IB[l,0] ,PD[l,0] ,RP[l,l] ,5[l,0] ,SP[l,0] ; 

*  42. 00  END; 


Figure  A4.  Network  Analysis  Program 
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MATCK  DATA  OUTPUT 


I 

Z 

IB 

PD 

RP 

S 

SP 

BRANCH 

OHMS 

AMPS 

WATTS 
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PERCENT 

0 
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0.08 

1.00 

0.08 

7.53 

1 

2.00 

1.U8 

U.38 

1.00 

4.38 

1+38.17 
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2.00 

-it.  60 

1+2.37 

1.00 
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1.00 
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5.00 
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1.00 
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1+.00 

1.91 

14.61+ 

1.00 

Ik. 6k 

11+61+.37 

Figure  A5»      Output  of  Network  Analysis  Test  Program 
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APPENDIX  B 


REGRESSION  MODEL  SUBROUTINE 


I   SUBROUTINE  ITERATE  LISTING 

COMMON  CCSAVIfc,20G) ,IDW0R0(74) 

DOUBLE  PRECISIUN  ZPRL150.50)  GNL0O030 
COMMON  NMAX,NNODE  ,  NTE RMS , NUMBL , NUMBR , NUMBC ,  IR TN , NTR AC E  ,  NSWTCH , KTG, GNLG005  0 

1  NPRINT(IO)  GNLG0060 
COMMON  E(2CC),EMIN(200),EMAX(2CC)  ,  AMP  (  200  )  ,  AMPM  I  N  (  200  )  t  AMP  MAX  (  200  )  GNL  000  7  n 

COMMON  YI 200),YMIN<2  00) , YM AX< 200 ) , N I N I T ( 200 ) , NF IN ( 200 ) , MODE  1 ( 200 )  GNL  000  tK) 

COMMON  Y TERM (200) , YTERMH<2  00) iYTERML(200I  ,  I ROWT ( 200 ) , I  COL T < 200 )  GNL00090 

COMMON  ERR0R1,ISEQ,MSEQ,M0,NUMM0, VFIRSTI50) , VSECNDI 50  I , VL AST ( 50 )  GNL 00 100 

COMMON  MOBRN(50),MOPARM(50) ,MOSTEP( 50),INX0UT(4)  GNL00110 

C  GNL00120 

C      THE  FOLLOWING  VARIABLES  ARE  USED  ONLY  IN  THE  ECAP  D.C.  ANALYSIS  GNL00130 

COMMON  AX1.SMLEPI 50) ,CURR ( 200 ), SMLE ( 200 ), EQUCUR ( 50 ) ,EX(200)  GNL00  150 

COMMON  EBI2CO),AMPX(200),AMPB<200),VNOM«50)  ,  STDSQ ( 5  0 )  ,L  ,  M  ,  IITOL  GNL 00  160 

COMMON  JX1,  JX4,JX5, DELTA, DUNK28)  GNL00170 

COMMON  NW0RCS<72) ,NMCO( 2 , 201 , KLABEL ( 41 , KPUNC(5) ,  INDCI2,20)  GNL001S0 

COMMON  INPUTEM9) ,NBC0(20) ,  KTYPE(5) , NBL ANK , NOEXEC , I TOL , NEQU I M , IPC  GNL 00 191 

COMMON  INVAL,LL,IC0L,LTYPE,KCCL,NCUIT,ITRANS,K0,KS,KELAST,NUM,M1  GNL 00 192 

COMMON  M2tM3,KCARD,KG,NP,NTR,MAC,HNODE,TNUM,NOEL,N0E,N0I,NOIC  GNL 00 193 

COMMON  EQUIVN(20J ,K0UT(2,10)  GNL00194 

COMMON  MATA(200,4,3),YX(200I,V8(200) , YTERMXI 200 ) , YT ERMB ( 200 )  GNL00210 

COMMON  WCMAX(50),WCMIN(50)  GNL00220 
COMMON   CCSAV(2,200) 

OOUBLE  PRECISION  SMLEP ,CURR, SMLE , EQUCUR  GNL00240 

C  GNL00270 
C      THE  FOLLOWING  VARIABLES  ARE  USED  ONLY  IN  SUBROUTINE  ITRATE 

DIMENSION  NITER(6,2),C0EFS150,5), ICRESI50) , IPTYPE ( 50) , I PLCC (50)  GNL 00280 

DIMENSION  ICURR<2  5),T0L(25>,AMP0LD(25),BETA(200),GM(200)  GNL 00290 

DOUBLE  PRECISION  MACURI200)  WNL00295 

REAL*8   IREF  WNL00297 
C 

C       INITIALIZE  AND  DEFINE  DICTIONARY  AND  INDICATORS  GNL00310 

IF  (NTRACE)  1,3,1  GNL00330 

1  WRITE  (6,2)  GNL00340 

2  FORMAT  (•  SUBROUTINE  ITRATE,  REGRESSION  MODEL  VERSION  ENTERED') 

3  NITERU.ll*  -482328512  GNL00360 
NITER(l,2)=  -985644992  GNL00370 
NITERI2,l)=  -1019199424  GNL00380 
NITER(2,2)=  -465551296  GNL00390 
NITER! 3,  I  )=  -482328512  GNL00400 
NITER13.2)*  -700432320  GNL00410 
NITER(4,1)=  -985644992  GNL00420 
NITERI4,2I»  -415219648  GNL00430 
NITER(5,1)«   -733986752 
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NITER<5,2)=       -700432320 
NITER(6,l»=      -650100672 

NITER<6,2)=      -985644992  GNL00440 

DO    300    1=  1. 5C  GNL00450 

ICR€SII)=0  GNL00460 

IPT^PEl 1 1=0  GNL00470 

IPLOC(l)=0  GNL00480 

DO    300    J=l»5  GNL00490 

300  COEf S(IiJ)=0.0  GNL00500 
DO  301  1=  It  25  GNL00510 
ICURRU>  =  0  GNL00520 
TOL(l)=0.C  GNL00530 

301  AMPOLD(I)=0.0  WNL0053? 
DO  102  1=1,200  WNL00534 
BETA(I)=0.0  WNL00536 
MACURII  )    =    0.0  WNL00538 

302  GM(I)=0.0  GNL00540 
INDCUR=0  GNLOO550 
INDTQL=0  GNL00560 
IN0TEM=0  GNL00570 
TEMP=0.0  GNL00580 
TEMPI=0.0  GNL00590 
TEMPF=0.0  GNL00600 
TEMPC=0.0  GNL00610 
IINMM=0  GNL00620 
ITEHP2=i  GNL00630 

C  GNL00640 

C               READ    AND    LIST    DATA    CARD  GNL00660 

1000  READ    15,4)    NWCRDS  GNL00670 

4  FORMAT  (72A1)  GNLOO6R0 
URITE    (6,5)     NWORDS  GNL00690 

5  FORMAT  (IX.72AI)  GNLOOfOO 
KCARD=KCARD*l  GNL00710 
NOEXEC=N0EXEC*NQUIT 

C 

C  CHECK    FOR    ■  IC    CARD  t#ww 

IF    (NW0RDS(l)-INDC(l,15))     1004,1001,1004 

1001  IF    CNtaORDSm-NMCDCltin     1 004,  1C02,  1004 

1002  DO  1003  ICOL=l,72 

1003  IDWORDI  ICOL)  =  NWORCS<  ICOL) 
GO  TO  100C 

C 

1004  CONTINUE                     _  ,An„  GNL00720 
IF  INWOROS(1)-INPUTBI8))  2 000 , 1 OOC, 2000  GNL00730 

C                                          ,  GNL00  740 

C       SUPPRESS  BLANKS  FROM  COLUMNS  7-72  GNL00750 

C  GNL00760 

2000    KCCL=6  GNL00770 

DO    7    IC0L=7,72  GNL00780 

IF    (NWORDSUCOL)-NBLANK)     6,7,6  GNL00790 

6  KCOt=KCCL»l  GNLOO^OO 
NWORDS(KCOL)=NWORDS(  ICOL)  GNL00810 

7  CONTINUE  GNL00H20 
C  GNL00H10 
C  CHECK  CONTENTS  OF  COLUMNS  1-6  GNL00840 
C  GNL00850 

CO    II     1C0L=1,5  GNLG0860 

IF    (NWOROSI  ICOD-NBLANK)    8,11,8  GNL00870 

8  DO    9    LTYPE=1,2 
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IF  INWORDS(  ICOL)-KLABEL(LTYPE))  9,200,9  GNL00380 

9  CONTINUE  GNL00890 

M3«5  GNL00900 

15  ITRANS-6  GNL00910 

10  WRITE  (6,3006)  Mi  GNL00920 
WRITE  (6,3007)  KCARD.ICOL  GNL00921 
WRITE  (6,3008)  GNL00922 
NQliIT=l  GNL00923 

3006  FORHAT  I'C  ****♦  ERRCR  NO.', 13,'  ****  MS6  FROM  SUBROUTINE  ITRATE') 

3007  FORMAT  (•    CARO  NG.=',I3,3K,   'APPROXIMATE  COLUMN  NO.  IS', 13, 
l'HSG  FROM  SUBROLTINE  ITRATE*) 

3008  FORMAT  (//)  GNL00926 
GO  TO  1000  GNL00930 

16  FORMAT  ( '0' , 'SUBROUTINE  ITRATE  ERRO«  MSG  -  CHECK  YCUR  INPUT  DATA  0 
1VER  CAREFULLY,  PGM  SAYS  YOU  HAVE  MADE  A  GOOF') 

11  CONTINUE  GNL00950 
IF  (KCOL-6)  500,12,13  GNL00960 

12  M3=4  GNL00970 
GO  TO  15  GNL0098D 

500  ITRANS=5  GNL00990 

GO  TO  10  GNL01JOO 

C  GNL01010 
C      CHECK  FOR  TEMPERATURE,  CURRENT,  TOLERANCE,  EXECUTE,  MCOIFY  OR 
C      RETURN  STATEMENT,   INDEX  NUMBERS  1,  2,  3,  4 ,  5,  6  RESPECTIVELY. 

C  GNL01030 

13  ICGL=7  GNL01040 
DO  20  IDhNT=l,6 

IF  <NWORDS(  ICCL)-NITER(IDENT,1) )  20,14,20  GNL0D60 

14  IF  INWORDSI  IC0L*1)-NITER1 IDENT,2))  20,25,20  GNL01T70 

20  CONTINUE  GNL01080 
M3=12  GNL01090 
GO  TO  15  GNLOUJO 

21  M3  =  34  GNLOUIO 
GO  TO  15  GNL01120 

25  GO  TO  ( 26,26,26, 800C, 4000, 9999) ,IDENT 

C  GNL01L40 

C      GET  PAST  •=•  SIGN  GN  TOLERANCE,  TEMPERATURE  OR  CURRENTS  GNL01LS~> 

C  GNL01160 

26  IC0L=ICCL*1  GNL01170 
IF  IICOL-KCCL)  27,21,21  GNL011R0 

27  IF  INWORDSI  ICCL  )-KPUNC( 1) )  26,28,26  GNL01190 

28  ICOk.=  lCOL*l  GNL01200 
GO  TO  (50,29,39),  IDENT  GNL01210 

C  GNL01220 

C      PROCESS  TOLERANCES  AND/OR  CURRENTS  GNL01230 

C  GNL01240 

29  IF  IINDCUR)  21,30,21  GNLOl25n 

30  IN0CUR=1  GNLC1260 
41  1=0  GNLG1270 

31  CALL  ECB09  GNL01280 
I=U1  GNL01290 
NUM^TNUM  GNL0  1300 
IF  (1-25)  42,42,21  GNL01310 

32  ICURR(I)=NUM  GNLG1320 
I T=I  GNL01330 

33  IC0L=IC0L»1  GNL01340 
IF  IICOL-KCCL)  34,34,1000  GNL01350 

34  IF  INWORDSt  ICGL  )-KPUNC(3) )  36,35,36  GNL01360 

35  IC0L=IC0L*1  GNL01370 
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GO    TO    31  GNL01380 

36    M3=34  GNL01390 

GO    TO    15  GNL01400 

39  IF    (INDTOL)     21, 40, 21  GNL01410 

40  INDT0L=1  GNL01420 
GO    TO    41  GNL01430 

42  GO    TO    (500,32,43),IDENT  GNL01440 

43  TOLU)  =  TNUM  GNL01450 
JT=I  GNL01460 
GO    TO    33  GNL01470 

C  GNL01490 

C  PROCESS    TEMPERATURE  GNL01490 

C  GNL01500 

50  IF    4IC0L-KCCL)     51,51,36  GNL01510 

51  IF    UNDTEM)    36,52,36  GNL01520 

52  INDTEM=l  GNL01530 
CALL  ECB09  GNL01540 
TEHPI=TNUM  GNL01550 
IC0L=IC0L*1  GNL01560 
IF    (IC0L-KC0L1     54,54,53  GNL01570 

53  TEMP=TEMPI*273.0  GNL01580 
ITEMP2=1  GNL01590 
GO    TO    1000  GNL01600 

54  IF    (NWOROSUCOD-KPUNCm  )     36,55,36  GNL01610 

55  IC0L=IC0L»1  GNL01620 
CALL  ECB09  GNL01630 
TEMPC=TNUM  GNL01640 
ITEMP2=2  GNL01650 
ICOL-ICOLU  GNL01660 
IF   .1  ICOL-KCCL)     56,56,36  GNL01670 

56  IF    INWORDSC ICOL )-KPUNC (4 I)     36,57,36  GNL01680 

57  IC0L=IC0L»1  GNL01690 
CALL  ECB09  GNL01700 
TEMPF=TNUM  GNL01710 
GO    TO    1000  GNL01720 

C  GNL01730 

C  GET    COEFFICIENTS    AND    OATA    TYPE  GNL0W40 

C  GNL01750 

200  IC0L=ICCL*1  GNL01760 
IF    IICOL-KCOL)    201,201,21  GNL01770 

201  I  IN4JM=  I  INUMM  GNL01780 
IF    IIINUM-50)    203,203,202  GNL01790 

202  M3=36  GNL01800 
GO    TO    15  GNL01810 

203  CALL  ECB09  GNL01820 
NUM*TNUM  GNL01830 
IPLOCd  INUM*=NUM  GNL01840 
IF    INWORDS<7)-INPUTB(5) I     204,207,204  GNL01850 

204  IF    CNWOKDS(7)-INPUTB(2) J     205,208,205  GNL01860 

205  IF    INWORDSI 7 )- INPUTBI 4) )    206,209,206  GNL01870 

206  IF    INWORDSm-INPUTBOJ)    225,224,225  GNL01880 

207  MTYPE=1  GNL01890 
GO    TO    210  GNL01900 

208  MTYPE=2  GNL01910 
GO    TO    210  GNL01920 

209  MTYPE=3  GNLC1930 

210  IPTYPEI IINUM)=MTYPE  GNL01940 

211  IC0L=IC0L*1  GNL01950 
IF    IICOL-KCCL)     212,212,21  GNL01960 
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212  IF  (NWORDSl  ICCL J-KPUNCI2) »  211,213,211  GNL01970 

213  ICOL-ICOLU  GNL01980 
IF  UCOL-KCOL)  214,214,21  GNL01990 

214  CALL  ECB09  GNL02000 
NUM-TNUM  GNL02010 
ICRESl IINUM)=NUM  GNLO2020 
ICOL=ICGL»l  GNL02030 
IF  4ICOL-KCCL)  215,215,21  GNL02040 

215  IF  (NWORDS<  ICOL  I-KPUNCI4II  206,216,206  GM02050 

216  ICCfc=ICOL*l  GNL02060 
IF  (ICOL-KCCL)  217,217,21  GNL02070 

217  IF  (NWORDSl  ICOL )-KPUNC ( l))  206,218,206  GNL02J80 

218  IICOL=0  GNL02090 

222  IC0L«IC0L*1  GNL02100 
IF  4  ICOL-KCCL  I  220,220,219  GNL02110 

219  IF  IIICOL)  5C0, 206, 1000  GNL02120 

220  IICCL=I1CCL*1  GNL02130 
IF  <IIC0L-5)  222,223,206  GNL02140 

223  CALL  ECB09  GNL02150 
COE€S( I  INUM.I  ICOL)  =  TNUM  GNL02160 
ICCL=IC0L*1  GNL02170 
IF  (ICOL-KCCL)  221,221,1000  GNL021B0 

221  IF  (NW0R0S(  ICOL  )-KPUNC(3) )  206,222,206  GNL02190 

224  MTY*E=4  GNL02200 
GO  10  210  GNL02210 

225  M3=36  GNL02220 
GO  TO  15  GNL02230 

C  GNL02240 

C      BEGIN  EXECUTION  OF  ITERATION  TECHNIQUE  GNL02?U>0 

C  GNL02260 

8000  IF  4NTRACE)  890  1 , 89C0, 890 1  GNL02270 

8901  WRITE  (6,89081  GNL02280 
WRITE  (6,8902)  ( TCL ( I  )  , I  =  1 , I T  )  GNL02290 
WRITE  (6,89031  (  ICURR ( I ) , I  =  1 ,  I T )  GNL02300 
WRITE  (6,89C6)  GNL02310 
WRITE  (6,8904)  ( ( COEF S ( I , J ) , J= I , 5 ) , I  =  1 , I  I NUM )  GNL02320 
WRITE  (6,8907)  GNL02330 
WRITE  (6,8905)  I  I  ,  I  CRES<  I  )  ,  I  P  T  YPE  (  I J  ,  IPLOC  (  I )  ,  I  =1  ,  20  )  GNL023<»0 
WRITE  (6,8911)  GNL02350 
WRITE  (6,8910)  (CLRR (J )  ,  J  =  1  ,NMAX )  GNL02360 

8902  FORMAT  I  IHO , 1 CF 1 0 .4 )  GNL02370 

8903  FORMAT  (  1HC , 1C ( 5X , I  5 ) )  GNL02380 

8904  FORMAT  (1X,5E17.7)  GNL02190 

8905  FORMAT  (  IX ,  I  3  ,  5X ,  I  5  ,  5X , I  5 , 5X , I  5 )  GNL02400 

8906  FORMAT  (  IHO , 1 2HC0EFF IC IE  NT S/ )  GNL02410 

8907  FORMAT  I  1  l-C  ,  2  X  ,  1H  I  ,  7X  ,  5H  IC  RES  ,4  X  ,  6HI  PT  YPE  ,  5  X,  5H  I  PL  OC/  )  GNL02420 

8908  FORMAT  (  1  HO  ,  2 3HT0LERANCES  AND  CURRENTS!  GNL02430 

8910  FORMAT  (1X,8E16.7)  GNL02450 

8911  FORMAT  (  1  HO  ,  81-CURR  ENT  S  /  )  GNL02460 
8900  IF  4  INDCUR+INDT0L-2)  80C1 , 8C02 , 8001  GNL024H0 

8001  M3=35  GNL02490 
GO  JO  15  GNL02500 

8002  IF  41T-JT)  8CC3, 8004, 8003  GNL02510 

8003  M3=36  GNL02520 
GO  TO  15  GNL02S30 

8004  CONTINUE  GNL02540 
NUM=0  GNL025bO 
GO  TO  (8098.8C5C) ,ITEMP2  GNL02560 

C  GNL0257O 


SINGLE  TEMPERATURE 

8008  CONTINUE 

60  TO  1000 

RANGE  OF  TEMPERATURES 

8050  TEMPC=(TEMPF-TEMPI )/TEMPC 
TEMPF=TEMPF*273.0 
TEMPI=TEMPI»273.0 
TEMP=TEMPI 
GO  TO  8098 

8053  TEMP=TEMP»TEMPC 
NUM=0 

IF  ITEMP-TEMPF)  8098, 8098, 8054 

8054  CONTINUE 
GO  TO  100C 


;      STORE  PRESENT  CURRENT  VALUES 

8098  DO  8099  1  =  1, IT 
J=ICURR<I  ) 

8099  AMPOLD(I)=CURR( J) 
NUM=NUM*l 

CALCULATE  NEW  DEPENOENT  VALUES 

8100  00  8105  1=1,  IINUM 
J=ICRES(I  ) 

IREF  =  10CCCCC.C 

CUTOFF  =  1.0 

ZOFf  =  5555555. CCCO 

MACUR(J)  =  <CURR( J) )*IREF 

SAVE    =    SNGLIMACURU) ) 

MTYPEMPTYPEd) 

K=lPLOC(l) 

NFBft    =     ICGLTIKI 


C 

C  C 
8101 
8920 


614 

616 

C  616 

8928 


GOTO    18101,8102,8103,8104) ,MTYPE 

ALCULATE    NEW    VBE 

IFISAVE    .LE.    CUTOFF)     WRl TE ( 6, 8920 )     J, SAVE, J 

FORMAT!'     «,20X, 'CALCULATION    OF    VBE',  5X    ,•    CURR(',I3,'I    =• 

1    F20.8,     5X,«FRCM    BRANCH    =',I3) 

IFISAVE     .LE.    CUTUFFI    GOTO    614 

ALOJR    =    ALOG(  SAVE) 

VALUE=COEFSI  I,  1)    ♦    ALCUR*(COEF S( I , 21 ♦    COEFS < I , 3 ) *ALCUR > 

EIKI       =    VALUE       ♦    CCSAVI3,K) 

GOTO    6  16 

VALUE  =  COEFSI  1,1) 

E(K1   =  VALUE 

CONTINUE 

WRITE(6,8928)J,SAVE,J,CURR< J)  ,  MTYPE , K ,NFBR , ALCUR , VALUE 

FORMAT!  'O'.IOX,  'CALCULATION  OF  VBE  •  ,  5X  ,  •  J=  •  I  3,  5X,  •  IC  =',F8.3, 
15X,'CURRN('I3,' )  =  'F12.6.5X, 'MTYPE  ='I3,5X,'K  ='13,/ 
2'  ',34X,'KFBR  =  •  I  3  ,  5X , 'LN( IC )  ='F12.6,5X,'  VBE  ='F12.6/  ) 

GOTO    8105 
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GNL02580 
GNL02590 
GNL02600 

GNL02620 
GNL02630 
GNL02640 
GNL02650 
GNL02660 
GNL02670 
GNL02630 
GNL02690 
GNL02700 
GNLO2710 
GNL02720 
GNL02730 

GNL02750 
GNL02760 
GNL02770 
GNL02780 
GNL02790 
GNL02800 
GNL02810 
GNL02820 
GNLC2830 

WNL02850 
WNL02851 


WNL 

WNL02853 

WNL02854 

WNL02855 

WNL02356 

WNL02860 

WNL02861 

WNL028X 


WNL02870 

WNL02872 

WNL0287 

WNL02873 

WNL02875 

WNL02876 


WNL02878 
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C   CALCULATE  NEH  BETA 
C 

8102  IFI6AVE  .IE.  CUTOFF)  WRI TE I  6, 8922  I  J, SAVE, J  MNL02880 
8922  FORMAT!'  ',  20X,  "CALCULATION  OF  BETA*,       5X  ,•  CURR(M3)'I  «•  , 

1  F20.8,  5X,«FRCH  BRANCH  »',<3I 
IF(SAVE  .LE.  CUTOFF)  GOTO  624 

ALCUR  ■  ALOSISAVE)  WNL02883 

VALUE  =  COEFSd.l)  ♦  ALCUR*«COEFS( I , 2 )  ♦  COEFS ( I , 3 ) *ALCUR)         WNL02B84 
Y(MFBR)  «  l.C/(ZOFF)  WNL02R87 

BETA(K)  ■  VALUE  WNL0288 

GOTO  626  WNL02885 

624  VALUE  =  0.0  WNL028R6 

BETA(K)  =  VALUE  WNL02388 

GOTO  626  WNL02889 

626  YTERM(K)  =  VALUE  *  Y(NFBR)  WNL02892 

C      WRITE (6,  8 929)  J, SA VE , J ,CURR ( J) , MTYPE , K ,NFBR , ALCUR , VALUE , K ,YTERM ( K I 
8929  FORMAT! '0' ,10X, 'CALCULATION  OF  BETA • , 5X, • J= • 15 , 5X , •  IC  =  'F8.3, 
15X, 'CURRNt ' 13,' )=  'F12.6.5X, 'MTYPE  ='I3,5X,'K  ='13,/ 
2'  «,34X,'NFBR  =  M3,5X, 'LN<  IC)  = 'F 12 .6,5X,  •  BETA  ='F12.6/ 
3'  ',33X,'YTERM(  •  12,  •  )  =',F12.6/) 
GOTO  8105  WNL02896 

C 
C 

C   CALCULATE  NEM  VALUE  OF  GM 
C 

8103  IFISAVE  .LE.  CUTOFF)  WRI TE (6, 8924 )  J, SAVE, J 

8924  FORMAT!'  •  ,  1  5X , • C ALCULAT I  ON  OF  GM   ',/'0',20X  , • ICURR ( ' , I  3  ,  • )  =• 
1  F12.6, 10X,'FR0M  BRANCH  =',I3) 
VALUE=C0EFS<1, 1)  ♦  ALCUR*  I  COEF  S  1 1  ,  2  )  ♦  COEFS  II  ,3  )*ALCUR  )  WNL03012 

GM(  K)=VALUE  WNL03014 

YTERM<  K)=VALUE  WNL03016 

GO  TO  8105  GNL03030 

C 

C 

C   CALCULATE  MEW  VALUE  CF   Y 

C 

8104  IFISAVE  .LE.  CUTOFf-  )  WRI  TE  I  6  ,  8926  )  J, SAVE, J 

8926  FORMAT!'  •  ,  1  5X , ' CALCULAT I  ON  OF   Y   ',/'0',20X  , • I CURR ( '  ,  I  3  ,  '  )  =' 
1  F12.6,  10X,  'FRCM  BRANCH  =',I3) 
VALUE=COEFS(  I,  1)  *■    ALCUR*  I  COEF  SII  ,  2  )  ♦  COEFS  <  I  ,3  )  *ALCUR  )  WNL03042 

V(NFBR)=1.0/VALUE  WNL03044 

8105  CONTINUE  GNL03050 
C 

C 

C      SOLVE  DC  PROELEP  GNL03070 

C  GNL03080 

DO    8106     I=*1,NMAX  GNL03110 

YX<I)=Y(I)  GNL03120 

EXU)=EII)  GNL03130 

8106  AMPXI I)=AMPII)  GNL03140 
IF  INTERMS)  ei07,81C9,8107  GNL03150 

8107  DO  8108  I=l,NTERMS  GNL03160 

8108  YTERMX!  I  )  =  YTERM( I  )  GNLC3170 

8109  CALL  ECB22I2PRL)  GNL031B0 
CALL  ECB23  GNL03190 
CALL  ECB26IZPRL)  GNL03200 
CALL  ECB25  GNL03210 
IF  (NUM-20)  £200,6300,8200  GNL03230 

C  GNL03240 
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CHECK  TOLERANCES 

8200  00  8201  I-ltIT 
J-ICURR(I) 
SAVE=CURR( Jl 

IF  (ABSl AMPQLD(I)-SAVE) 

8201  CONTINUE 

;      OUTPUT  RESULTS 


-  A8S(AMP0LD(I)*T0LII)))  820  1 ,820  1 1  8098 


8302  WRITE  (6,8302)  NUM 

8303  FORMAT  (•-•,  'NUMBER  OF  ITERATIONS  =  • ,  1 3/1 
IF  (INDTEM)  8202,8203,8202 

8202  TEMPA=TEMP-273.0 
WRITE  (6,9000)  TEMPA 

9000  FORMAT  ( 1H0 , 14HTEMPERATURE  *  ,F10.3/» 

8203  00  8204  I=1,NMAX 
YX(I)=Y( I ) 

EX(I )=E< I ) 

8204  AMPKd  )=APP(  I) 

IF  (NTERMS)  e205,82C7,8205 
820b  DO  8206  1=1, NTERMS 

8206  YTERMXl I)=YTERM( I ) 

8207  CALL  ECB22(ZPRL) 
CALL  ECB23 

C    THE  FOLLOWING  CODE  RECOGNIZES  AND  EXECUTES  THE  PRINT  NV,MI  CONTROL 
C    INPUT  DATA  STATEMENT.  RESULTANT  OUTPUT  IS  THE  DC  AND  NL  MATRICES. 

IFINPRINT(IO) )  107,107,108 
108  JX1=1 

CALL  ECB24IZPRL) 
107  CALL  ECB26WPRL) 

NTR=1 

IFINPRINT(IO)     103,103,104 
104    JX1=2 

CALL    ECB24UPRL) 
103    CALL    ECB25 

NTR=4 

WRITE    (6,9001) 

WRITE  (6,9002)  (  I  ,  ICCLT( I ) , IROWT( I ) ,BETA( I ) ,GM(  I ) , YTERM  (  I )  , 
1  1=1, NTERMS) 

9001  FORMAT  (//,IX,9H  T-NUMBER , 7X ,6H  IC0LT,11X,6H  IROW T, 6X ,4HBET A, 1 IX, 
12HGM,llX,5HYTERM/) 

9002  FORMAT  (  3X  ,  I  4  ,  1 IX  ,  I  4  ,  13X  ,  I  4,2F  I  5.  7,  E  1  5.  7) 
GO  TO  (8008,8053), ITEMP2 

8300  WRITE  (6,63C1) 

8301  FORMAT  (IHC34HN0  CONVERGENCE  AFTER  20  ITERATIONS/) 
GO  TO  8  30  2 

C 

C      BEGIN  MODIFY  ROUTINE 

C 


4000  CONTINUE 

GO  TO  1000 


9999 


EXIT  TO  MAIN  ROUTINE 

RETURN 
END 


GNL03250 
GNL03260 
GNL03270 
6NL03280 


GNL03300 
GNL03320 
GNL03330 
GNL03340 
GNL03350 
WNL03360 
GNL03370 
GNL03380 
GNL03390 
GNL03400 
GNL03410 
GNL03420 
GNL03430 
GNL03440 
GNL03450 
GNL03460 
GNL03470 
GNL03480 
GNL03490 


GNL03530 
GNL03540 
GNL03550 
GNL03580 
GNL03b90 
GNL03600 
GNL03630 
GNL03640 
GNL03650 
GNL03660 


GNL03670 
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2   COMMENTS  ON  PROGRAM  CHANGES 

C  I.  ORIGINAL  VERSION  OF  THIS  SUBROUTINE  MAS  WRITTEN  BY  G.  GUTH. 

C  CHANGES  MADE  BY  0.  WOOO  IN  INPUT-OUTPUT  ROUTINE  TO  PERMIT 

C  PARAMETER  VARIATION  ANALYSISt  USE  OF  VARYING   REFERENCE 

C  CURRENTS,  IREF,  AND  MODIFICATION   OF  THE  BASIC  ROUTINE  TO 

C  PERMIT  IT  TO  HANDLE  BOTH  CUTOFF  AND  ACTIVE  TRANSISTOR 

C  OPERTAING  REblONS. 

C 

C  2.  INSERTED  MACUR(J)  TO  IMPROVE  THE  FLEXIBILTY  OF  REGRESSION 

C  MCDEL.   THUS  BY  THE  USE  OF  A  SCALE  FACTOR  OF  10  =  IREF,  THE 

C  CURRENT  NORMALIZATION  INTHE  REGRESSION  MODEL  MAY  BE  CHANGED 

C  BY  A  FUTURE  MODIFICATION  PERMITTING  IREF  TO  BE  READ  IN  AT 

C  EXECUTION  TIME. 

C 

C  3.  CUTOFF  IS  THE  COLLECTOR  CURRENT  VALUE,  IC ,   AT  WHICH  BETA  IS 

C  SET  EQUAL  TO  ZERO.  CUTOFF  IC  IS  THAT  VALUE  OF  IC,  DETERMINED 

C  FROM  THE  REGRESSION  MCDEL  EQUATIONS,  FOR  WHICH  BETA  =  0. 

C 

C  4.  ZOFF  IS  ThE  IMPECANCE  LEVEL  SWITCHED  INTO  THE  BASE  BRANCH 

C  WHENEVER  THE  COLLECTOR  CURRENT  BECOMES  .LE.  CUTOFF. 

C  ZOFF  SET  EQUAL  TO  =  5555555. C  IN  PROGRAM,  THUS  ENABLING  THE 

C  USER  TO  SPOT  IN  THE  IMPEDANCE  LISTINGS,  EXACTLY  WHEN  THE 

C  PROGRAM  PRECICTED  ENTRY  INTO  THE  CUTOFF  MODE  OF  OPERATION. 

C 

C  5.  DUE  TO  THE  THE  UNIT  OF  CURRENT,  AMPERES,  ASSUMED  BY  THE  ECAP 

C  SUBROUTINES,  SET  IREF-1.0,  WHEN  A  1.0  AMPERE  REFERENCE  IS 

C  DESIRED.  IF  MICROAMPERES  IS  DESIREO,  SET  IR EF  =  1000000.0,  AND 

C  IF  MILLIAMPERES  IS  DESIRED,  SET  IREF=1000.0. 

C 

C  6.  THE  ARRAY  IN  COMMON,  CDSAV<2,200)  IS  USFD  TO  PASS  YTERMI200) 

C  AND   BETAI200)  VALUES  TO  NLAP  OUTPUT  ROUTINE  ECB25. 

C  CCSAV(1,J)  =  YTERM(J)  AND  CDSAVI2.K)  =  BETAIK). 

C 

C  7.  THE  FOLLOWING  INFORMATION  IS  TO  ASSIST  PGMR  IN  ADDING  MODIFY 

C  OPTION  TO  THIS  SUBROUTINE,  TO  ELIMINATE  PRESENT  REQUIREMENT 

C  FOR  DUPLICATING  CATA  DECKS  FOR  THE  CASE  OF  PARAMETER 

C  VARIATION  OPTION  RUNS* 

C  A.  MODIFY  ROUTINE  IS  ECA30 

C  B.  M  =  MODIFY  NO.,  NUMMO=  NO.  MOOS,  MOSTEP=  NO.  MOD  STEPS 

C  C.  MCPARM(M)  IS  KIND  OF  PARAMETER  TO  BE  MOHIFIEO 

C  D.  MOBRN(M)  IS  CARD  SEQUENCE  NUMBER 

C  E.  NTERM  IS  TOTAL  NUMBER  OF  DATA  CARDS 

C  F.  VFIRSTIM)   CONTAINS  VALUE  OF  ALL  PARMS  TC  BE  MODIFIED 

C 

C  8.  NOTE  THAT  GM  AND  Y  NEW  VALUE  COOES  HAVE  NOT   BEEN  MODIFIED. 

C 

C  9.  THIS  SUBROUTINE  AS  NOW  CODED  WILL  GIVE  ACCURATE  RESULTS  FOR 

C  THE  NL  OPTION  OF  NLAP,  PROVIDED  THAT  THE  USER  IS  VERY 

C  CAREFUL  TC  KEEP  INPUT  DATA  SO  THAT  THE  CIRCUIT  OOES  NOT  GO 

C  DEEP  INTO  THE  CUTOFF  REGION  OF  OPERATION,  OTHERWISE  CURRENT 

C  CONVERGENCE  MAY  NOT  BE  OBTAINED  DUE  TC   PRESENT  REGRESSION 

C  MODEL  EQUATION  FOR  VBE  NOT  HAVING  AN  EXTERMAl  CIRCUIT 

C  CURRENT  VARIABLE  INCLUDED  AS  AN  INDEPENDENT  VARIABLE. 
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APPENDIX  C 


LINEAR  REGRESSION  PROGRAM 


I   LINEAR  REGRESSION  PROGRAM  LISTING 


COMMON    A, ACHCDLD,  ICDCT 

COMMON  AI(40),ANAM(25,40),FMT<18) ,KZ( 20) , KCHG , T I TLE ( 16 ) 

DOUBLE  PRECISION  A  (  40  )  ,  AC  HG  (  4  C  ,  4 )  t  DLO  ,  I  COC  T 

DOUBLE  PRECISION         AI 0 ( 20  )  ,  CS ( 20 )  ,C ( 2C , 20 ) , V ( 20  I  ,  BB I  20 ) 

DOUBLE  PRECISION  TST,      SSS, STDV,SPST 

DOUBLE  PRECISION  AMEAN ,C0, AKNT , SOUT , SIGN ,PCT, STER , AST 

DOUBLE  PRECISION  BETA, BVAR ,BSTD, TRAT , SIZE, AM, DF ,RR , R , ADJR, ADJRR 

REAL    DSAVUC101) 


1056 

FORMAT  (18A4) 

1000 

F0RMATI6X.I2, 

510 

FORMAT  { 6  X  ,  24 

1001 

FORMAT  (6X,16 

1057 

FORMAT  <«-',2 

1002 

FORMAT   C»l»« 

8001 

F0RMATI6X.F2. 

8002 

FORMAT  ('C', 

615 

F0RMAT('-',2X 

9000 

FORMAT  HHOtl 

9002 

FORMATi 1H  ,7F 

1066 

FORMAT  <'-',2 

1 

L2X,    'SAMP  S 

1067 

FORMAT  <•  ',1 

1006 

FORMAT  (6X.10 

488 

FORMAT  (44HCT 

LLE  SIZE  =  ,/,l 

2000 

FORMAT  (8H0VA 

446 

FORMAT  I'C'/i 

459 

FORMAT  (•  ',3 

455 

FORMAT  <•  ',6 

4002 

FORMAT  I'2',I 

660 

FORMAT   CO', 

3009 

F0RMAT('-',18 

490 

FORMAT! '0'15X 

3010 

FORMAT  (  'C',6X 

IF8.3  ) 

472 

FORMAT( «0',2X 

LR  OF  EST  *•  , 

475 

FORMAT!'-',  2 

3012 

FORMATCO',  2 

L,3X,'T-RATI0' 

3053 

FORMAT  <•  ',2 

469 

FORMAT  (1HC,/ 

3091 

F0RMATC»-',2X 

IX, 12, IX, 3(11, IX), F9. 01 

Al) 
A4) 
X, 'INPUT  DATA  FORMAT  =  '.2X.16A4    > 

2X,  18A4) 
C,1X,F2.0,1X,F2.0,1X,F6.0  ) 

2X,  'ACHGJ ' ,12,'  )  =  «,F3.0,1X,F3.0,1X,F3.0,1X,F6.0   I 
, 'INPUT  OATA',   4X,15(9A1,3X)I 
CH  CJ  I,J)=    ,6F6.2,//) 

SIZE=   .F10.21 
X.'VAR  N0.',6X,'NA*E',17X,'MEAN',3X,   'POP  STD  DEV, 
TD  DEV   ) 
5.5X.2CAI,  3F12.5J 
(Fl.O,  1X1) 

00  FEW  DEGREES  OF  FREEOOM,  TRY  ANOTHER  EQN , /,  14H  SAMP 
8HINDEP  VAR  IN  EQN  «,F5.0) 

R  N0..I3.20HLEA0S  TO  S  INGULAR I  TV .5X8HVALUE  IS.F12.8) 
•0',2X,  'INTERCORRELATION   MATRIX'    ) 
SX^^./VO' ,24X,9I5) 
Xt20Al,5X,I4,2X, 16F6.2) 
1 • 1  IS!  •-  •  )  ) 

2X,18A4  ) 
X, 'DEPENDENT  VARIABLE  IS  NUMBER   -',13,'  -»,3X,  20A1) 

, 'SAMPLE  SIZE  =«,F6.1,7X, 'DEGREES  FREEDOM  ='F6.1   I 
,'R**2  =  •  ,F5.  3,7X, 'R.  =  • ,F5.3 , 8X, • STC  ERR  OF  EST  =•  , 

,'ADJ  R**2  =  • ,F5.3,3X, 'ADJ  R  =  • , F5 .3  ,4X, • AOJ  STD  ER 

F8.3     ) 

X, 'LINEAR  REGRESSION  EQUATION  COEFFICIENTS'   > 

X'VAR  N0.'4X'NAME'19X6HREG  WT,4X,'BTA  WT',4X,'WT  SIG' 

) 

X  ,  14, 3X,  21A1,  F12.5.3F11.5  ) 

//) 

,  'RESULTS  ARE  IN  RAW  SCCRES'I 
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3086  FORMAT  l'C',2)!,   'END  Of  PROBLEM'//  *2't 
430  FORMAT!'-',  2X,'END  OF  REGRESSION  PG«,,/'1,1 
TOL-0. 00005 
10  00  44  1=1,20 
Aid  I  )*0.C 
00  14  J«l,20 
14  C(  It J)=O.C 
C   REAO  TITLE 

REAO  ( 5,  1C56)  (TITLfcl  I  )  ,1=1,  16) 
C   REAO  NUMBER  OF  VARIABLES  CARD  ,  NVAR  ■  10  MAX 
16  REAO  (5,10CC  )  N  I  NO , N VAR , 1  COR , I  PR T ,  ICHG, T ST 
NIN0*NIN0*1 
ICOCT«0 
NVARM=NVAR*1 
M=NVARM 

1FINVARM  .LT.  1)  GOTO  650 
20  DO  24   1=2,  NVARM 
C   REAO  IN  NAMES  OF  DATA  VARIABLES  CARO  CCLS  7  -  27. 

24  REAO  ( 5,51C»( ANAMI J,I» ,J=1,24) 
C   READ  IN  DATA  FORMAT 

READ  (5, 10C1HFMTI  I  1,1  =  1,16  ) 
C   WRITE  TITLE  OF  OUTPUT  DATA 

WRITE  lb,   1CC2M  TITLE  U)  ,K=1,  18) 
C   WRITE  FORMAT  USED  TO  READ  IN  AND  WRITE  OUT  THE  PGM'S  INPUT  DATA 
WRITEI6.1C57) IFMT(K) ,K=1,16  ) 
IF!  ICHG  .LE.  C)  GOTO  3e 
KCH6=0 
NTRANS  =  4C 
DO  35    L  =  l,  MRANS 
C   READ  IN  DESIRED  VAR.  TRANSFORMATIONS,  MAX  =  10. 
REAO  (5,80C1  )  (ACHGIL.J), J=l,4) 
IFlACHG(L.l)  .LE.  0)  GOTO  3e 
C   WRITE  OUT  TRANSFORMATIONS 

WRITE  (6,8002)  L , ( ACHGI L , J )  ,  J=l ,4 ) 
KCHG=KCHG*1 
35  CONTINUE 
C   38  WRITE(6,8C03) 
C   WRITE  NAMES  CF  INPUT  DATA  VARIABLES. 

38  WRITEI6.615)  (( ANAMI JK , IJ )  ,  JK= 1 , 9)  , IJ=2 , NI NO ) 
C.  READ  IN  NIND    COLS  OF  DATA  PER  DATA  FORMAT  CARO,   REF.  STMT  1001 
40  REAO  I5,FMTM  A(  I)  ,  1  =  2, NIND) 
IFIAI2)  .GE.  TST)  GOTO  60 
C   WRITE  LIST  OF  INDEPENDENT  VARIABLE  INPUT  DATA. 
WRITE  16.FMT  )  ( A(  I  I  ,  I=2,NINDI 
A(  1)=1.0 
ICDCT=ICDCT*1 

IFUCHG  .GT.  C)  CALL  CHANGG 
DO  55  1=1, NVARM 
00  55  J=I .NVARM 
55  CI  I,J)=C(  I,JJ*A( I  )*AI J) 

GOTO  40 
60  SIZE=CI1,1) 
C 

C  WRITE  ROW  OF  ♦*  IF  A(2)  .GE.  TST,  PROVIDED   LAST  DATA  CARD  HAS  ALL  9'S 
C      OTHERWISE  THE  LAST  DATA  VARIABLE  IS  LISTED  IN  THE  OUTPUT  TWICE. 
WRITE  (6,FMT)(A( I ),I=2,NINDI 
Vt ll=OSQRT(C<  1,1)) 
DO  74   1=2, NVARM 
IF  IC(I.I)  .LE.  C)  GOTO   72 
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V(II«(SIZE*C<  I,  I  >-C(  1,  l)*CU,I>  )/<SlZE*SIZE) 

IF   IV(1)    .LE.    01    GOTO      72 

V1IJ=0SQRT(VU)>*V<1) 

GOTO    74 
72    ViIJ-l.0*vm 
74    CONTINUE 

SSS=0SQRT(SIZE) 

SSS=1.0/SSS 
WRITE    DATA    STATISTICS    HEADING 

WRITE<6, 1066) 

DO    76    I=2,NVASM 

STDV*VI I )*SSS 

SPST  =  V( I)»SSS 

SPST=SPST*SPST 

C0»ISIZE*SIZE)/((SIZE-1.0)*(  SIZE-  1.0)  ) 

SPST=SPST*CE 

SPST=DSQRT(6PST) 

AMEAN=CI  1* II/SIZE 

K=I-1 
WRITE  CATA  STATISTICS 

WRITE  (6,1C67)K,(ANAH( J, I ) ,  J  =  1  ,  Z 0 )  ,  AME AN , S TD V , SPST 
76  CONTINUE 

REAC  MODEL  DEFINATION  CARDS  (INVERT  CARDS) 
80  READ  I5,1CC6)(AH  I )  ,  l  =  2,NVARM ) 

All  11=1.0 

IF  (AH2I-3.C)  82,eCC,e50 
82  AKNT=0.0 

DO   84  I=1,NVARM 

IFU(I)  .NE.  1.0)  GOTO   84 

AKNT  =  AKM»1.C 
84  CONTINUE 

IF  (SIZE  .GT.  AKNT)  GOTO   e6 

WRITE  <6,48fi)SIZE,AKNT 
GOTO  80 
86  DO   92   I=1,NVARM 

IF    (AID(  I  )     .GT.    0)    GOTC    90 

csm=i.o/vm 

GOTO  92 
90  CS(I)=V( I) 
92  CONTINUE 

DO  94    I=1,NVARM 

DO    94  J=I,NVARM 
94  C(  I,J)  =  C(  I,J1*CS( I )*CS(J) 

DO  180  N=1,M 

IF  (AKN)  .EQ.  1)  GOTC  96 

SOUT=0.0 

GOTO  98 
96  S0UT=1.0 
98  IF  (SOUT  .EC.  AID(N))  GOTO  180 

IF(C(N,N)  .LE.  0)  GOTO  170 

IF  (CIN.M  .LE.  TCL)  GOTO  170 

101  C(N,N)=1.C/C(N,N) 

IF  ISOUT  .GT.  0)  GOTO  106 

102  DO  105  1  =  1, K 

IF  (I-NI103,1C5,1C4 

103  C(  I.NI=(  1.0-2.0*AID(I ))*C( I ,N) 
GO  TO  105 

104  C(N,I)=<1.C-2.0*AID( I))*C(N,I ) 
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105  CONTINUE  y 

106  00  111  I«1VM 

00  111  J=l,« 

SIGN*1.C-2.0»CABS< S0UT-AID(1)*AID( J) 1 

IF  -1  l-N  1107,111,110 

107  IF  IJ-NI108, 111,109 

108  C(  I, J)=C(  I,  J)*SIGN*C( I ,N)*C< J,NI*C(N,N) 
GO  10  111 

109  C(  I,  J)=C<  I, J  J ♦SIGN*CI I,N)*C(N,JI»CIN,N) 
GO  TO  111 

110  C(I,J)=C(I,J)*SIGN*C(N,I)»C(N,J)*C<N,N) 
HI  CONTINUE 

IF  (SOUT  .LE.  C)  GCTO  116 

112  DO  115  1=1, M 
SIGN=1.0-2.0*AIC( I) 
IF  I  I-N)113,115,114 

113  C(  I,N)=SIGN*C( I ,N)*C(N,N) 
GO  TO  115 

114  C(N,  I)  =  SIGMC(N,1  )*C(N,N) 

115  CONTINUE 
GO  TO   169 

116  00   168  1  =  1, M 

IF  (I-N)  161,  168,  162 

161  C(  I,N)  =  C(  I,M*C<N,N) 
GO  TO   168 

162  C(N,  1  )=C( N,  I)*C(N,N) 

168  CONTINUE 

169  AIC(N)=S01T 
GOTO  18C 

170  K  =  N-1 

WRITE  (6,200C)K,C(N,N) 

AI(N)=2.0 

GOTO  190 
180  CONTINUE 
190  DO  250   I=1,NVARM 

IF  (AICl  I  »)25C,241,242 

241  CS( I  J  =  V  I  I  ) 
GOTO  250 

242  CS(  I  )=1.0/V(I  ) 
250  CONTINUE 

IFllCOR  .LE.  01  GCTO  444 
WRITE  TITLE  •  INT ERCORRELATION  MATRIX* 
WRITE  (6,446) 
MKB=M-1 

WRITE  (6,459)1  IJ,  U=1,MKB) 

00  447  11=1, M 

00  447  JJ=II,M 
447  C( JJ,III=C( ( I.JJ) 

DO  <»50  11  =  2, P 

MK=1I-1 
WRITE  INTERCORRELATION  MATRIX  CATA. 

WRITE  ( 6,4  55  )(ANAM( JN.II) , JN=1,20) , MK  ,  ( C ( I  I  ,  J J ) ,JJ  =  2,II) 
450  CONTINUE 

00  458  11  =  1, N 

DO  458  JJ=1,M 
45  8  C(H,JJ)=C(  (I,JJ)*CS(  I  I)*CS(  JJ) 

ICOR  =  0 
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GOTO  80 
444  00  254   l»l,NVARM 

00  254   J-I.NVARM 
254  C(I,J)=C(  I,J)*CS< I )*CS( J) 

AN«0.0 

DO  299   I^ltNVARM 

IF  4AI01I  1-1.0)299,280,299 
280  AM=AH*1.0 

299  CONTINUE 

00  300   I*2,NVARM 

IF  4AKII  .LE.  1.0)  GOTO  300 
301  PCT  =  C( I,II/(V( I  )*V(  l>) 

STER-CI  I,  I  1/SIZE 

AST=SIZE/(SIZE-AM) 

AOSSR=STER*AST 

STER=DSQRT(STER) 

ADSSR=SQRT(ACSSR) 

OF=SIZE-AM 

RR=l.O-PCT 

R  =  DSQPT(fiH) 

ADJRR=1.0-(  l.C-RR)*( ( S I ZE-1.0 )/( SIZE-AM) ) 

IF  (AOJRR  .GT.  0.0)  GOTO  340 

ADJRR=0.0 

ADJR=0.0 

GOTO  344 
340  ADJR=DSQRT(ACJRR) 
344  K=I-1 
C   WRITE  REGRESSION  EQUATION  STATISTICS 

WRITE  <6,4002)IPRT 

WRITE  (6,  66C)  (TITLF(K),  K  =  1  ,16) 

WRITE  (6,3009)K,(ANAM(MK,l ) ,MK=1,20) 

WRITE  (6,490)SIZE,0F 

WRITE  (6,301G)KR,R,STER 

WRITE  (6.472JADJRR, AOJR.ADSSR 

WRITE(6,475) 
C   WRITE  TITLES  FOR  REGRESSION  EQUATION  COEFFICIENT  TABLE. 

WRITE  (6,3012) 

KJ=0 

00   320  J=1,NVARM 

IF    IAI0IJ4     .NE.    1.0)    GOTO    320 

IF   A  I-J1324,320,325 

324  BETA>Ct I.J) 
GOTO    326 

325  BETA=C« J, I) 

326  BVAR=C( I, I  )*C (J , J ) / < SI ZE-AMI 
BSTO  =  DSQRT(BVAR) 
TRAT=BETA/BSTD 

K=J-1 
KJ"KJ*1 
KZ(KJ)=K 
6B(KJ)=BETA 
BVAR=BETA*V( J)/V(I) 
C   WRITE  REGRESSION  EQUATION  COEFFICIENTS 

WRITE  (6,3053)K,( ANAH(NJ.J) ,N J=l , 21 ) , BETA ,BVAR , BSTO , TR AT 
320  CONTINUE 

300  CONTINUE 
C 

WRITE  (6,3091) 
GOTO  80 
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800    WRUE     (6,3066) 

IF(NVARM    .IE.    0)    GOTO    850 
820    GOTO    10 
850    WRITE(6,430) 

STOP 

END 

BLOCK       DATA 

COMMON  A, ACHG.CLD, ICDCT 

COMMON    AI440),ANAM(2  5,40),FNT(18)»KZ<20),KCHG,TITLE(18) 

DOUBLE  PRECISION  A( 40 I , ACHG44C ,*) ,0L0, ICOCT 

DATA  ANAM/'C',,0,,'N,,,S,,,I»  ,»A» f»M»,«T,,»  • , • T • , ' E« , »R' , • M' 
1  987*'  •/ 

END 

SUBROUTINE    CHANGG 

COMMON  A,ACHC,DLD,  ICOCT 

COMMON    A  I  140),  ANAfM  25,40)  ,FMT<  18)  ,K2(20)  ,  KCHCT  I  TLE  (  18  I 

DOUBLE    PRECISION    AI 401 , ACHGC4C,4> »DLD, ICDCT 

00    100    l=l,KCHG 
KBAB=ACHG(  1,1) 


N0=4DINT(  ACHC(I,2)  ) 

MD=IDINT(  AOCU.3)  ) 

LD=IDINT(  ACt-GII.A)  ) 

MD*MD*1 

N0=ND*1 

LD=LD*1 

DLC=AChG(  I  ,<.) 

GO  TO  <  1,2, 3, 4,5, 6, 7, 8, 9, 10, 11,12), KBAB 

1  A(MC)=A(N£)**CLO 

GO  TO  100 

2  A(MO)=1.0/(MND)**0LC) 

GO  TO  100 

3  A(MC)=A(N£)*A(LD) 

GO  TO  100 

4  A(MC)=AINC  l/MLC) 

GO  TO  100 

5  A(MO)=A(N0)-AILD) 

GO  TO  100 

6  A(MD)=DSQRTIA(ND) ) 

GO  TO  100 

7  A(MC)=DLOG(A(ND)) 

GO  TO  100 

8  AIMO)  =  0EXP(A(NC)  ) 

GO  TO  100 

9  A(MO)=A(ND)*A(LD) 

GO  TO  100 

10  A(MC)=ICDCT 

GO  TO  100 

11  A(MD)=A(NC)*OLD 

GOTO  100 

12  A(M0)=A(NC)*<DLD) 

100  CONTINUE 

RETURN 

END 
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2   TYPICAL  OUTPUT 

C3  HILLIAMPERE  LINEAR  REGRESSION  PGM,  C3  IC-NOHINAL  DATA,  VCE  -  2.0  VOLT 

INPUT  DATA  FORMAT  =    I16X,3F10*3I 
ACHGt  1)  =   7.   1.   «.     CO 
ACHGI  2)  =   3.   4.   5.      4. 


INPUT  CATJ 

i     COLLECTOR 

BETA 

VBE 

0.106 

0.0 

0.658 

C.540 

108.375 

0.702 

l.oeo 

108. OOC 

0.720 

1.293 

106.500 

0.  725 

1.616 

10  7.667 

0.730 

2.1*0 

104.800 

0.738 

2.658 

103. 6CC 

0.745 

3.163 

105. OOC 

0.749 

3.688 

101. OOC 

0.754 

A. 192 

100.800 

0.758 

5.213 

102. IOC 

0.763 

6.  176 

96.300 

0.770 

7.207 

10J. IOC 

0.772 

8.230 

1C2.30C 

0.775 

***♦*»♦*»»♦*»****#*♦♦*******»* 


INPUT  DATA  STATISTICS 


VAR 

NO. 

NAME 

MEAN 

POP 

STD  DEV 

SAMP  STD  D 

I 

COLLECTOR  CURRENT 

3.38018 

2.44  564 

2.63377 

2 

BETA 

96. 39586 

26.92444 

28.99555 

3 

VBE 

0.73993 

0.03071 

0.03307 

<. 

LN( IC,  MILLIAMPS) 

C. 79623 

1.12882 

1.21565 

5 

LN( IC,  MILLIAMPS)**2 

1.90821 

1.64586 

1.77246 

1NTERC0RRELAHCN   MATRIX 


COLLECTOR  CURRENT 

BETA 

VBE  3     0.86   C.67   1.00 

LNIIC,     MILLIAMPS)  4  0.86       0.68       1.00       1.00 

LNIIC,     MILLIAMPS)**2  5  C.59    -0.59       O.H       C.12       1.00 


1 

2 

1 

1.00 

2 

0.29 

l.CO 

3 

0.86 

C.67 

4 

0.86 

0.68 

5 

C.59 

-0.59 
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C3  MILUAMPERE  LINEAR  REGRESSION  PGM,  C3  IC-NCMINAL  OATA,  VCE  ■  2.0  VOLT 

DEPENDENT  VARIABLE  IS  NUMBER   -  18  -    BETA 

SAMPLE  SIZE  =   14.0        DEGREES  FREEDOM  =  11.0 

R**2  =  O.Sie        R  =  0.956         STD  ERR  OF  EST  =    7.705 

ACJ  R**2  =  C.9C3    ADJ  R  =  0.950     ADJ  STD  ERR  OF  EST  =    8.692 

LINEAR  REGRESSION  EQUATION  COEFFICIENTS 

VAR  NO.     NAME  REG  WT  BTA  WT  WT  SIG  T-RATIO 

0    CONSTANT  TERM  103.31701  3.83729     3.79670  27.21229 

4  LNUC,  MILLIAMPS)  18.08139  0.75807     2.07263  8.72391 

5  LNUC,  MILLIAMPS>**2  -11.17180  -0.68292     1.42151  -7.85908 


CEPENCENT  VARIABLE  IS  NUM8ER   -  18  -    VBE 

SAMPLE  SIZE  =   14.0  OEGREES  FREEDOM  =   11. 0 

R**2  =  0.999        R  =  l.OCO  STD  ERR  OF  EST  =    0.001 

ADJ  R**2  =  0.9S9    ADJ  R  =  l.OCO  ADJ  STD  ERR  OF  EST  =    0.001 

LINEAR  «EGRESSICN  ECUATION  COEFFICIENTS 

VAR  NO.     NAME  REG  MT     BTA  WT     WT  SIG    T-PATIO 

0    CONSTANT  TERM  0.71791    23.37929     0.00038  1879.57800 

4  LNUC,  MILLIAMPS)  0.02716     C. 99830     0.00021   130.24236 

5  LNUC,  MILLIAMPS)**2  0.00021     0.01114     0.00014     1.45294 

RESULTS  ARE  IN  RAM  SCORES 
END  OF  PROBLEM 
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C3  MICROAMPERE  LINEAR  REGRESSION  PGM,  C3  IC-NGMINAL  DATA,  VCE  *  2.0  VOLT 

INPUT  DATA  FORMAT  =    <16X,3F10»3) 
ACHG(  II  =   7.   1.   4.     C.O 
ACHGI  2)  =   3.   4.   5.      4. 


NPUT    CATA            CCLLECTCR 

BETA 

VBE 

106.500 

0.0 

0.658 

540. COO 

108.375 

C.702 

1080.000 

108. DOC 

0.720 

1292.  999 

106.500 

0.725 

1615.999 

107.667 

C.730 

2139.999 

IC4.800 

0.738 

2656.000 

1C3.60C 

0.745 

3183.000 

105. OOC 

0.749 

3688.000 

101. OOC 

0.754 

4191.996 

100.800 

0.758 

5212.996 

102.100 

0.763 

6175.996 

96.30C 

0.770 

7206.996 

103. IOC 

0.772 

8229.996 

102. 30C 

0.775 

****************************** 


INPUT  CATA  STATISTICS 


VAR 

NO. 

NAME 

1 
2 
3 
4 

COLLECTOR    CURRENT 

BETA 

VBE 

LN(  IC,    MICROAMPS) 

MEAN    POP  STO  OEV  SAMP  STD  DEV 

338C.  17693   2445.64334  2633.76975 

96.39586     26.92444  28.99555 

C. 73993      0.03071  0.03307 

7.70398      1.12882  1.21565 

LNIIC,  MICR0AMPS)**2     60.62560     15.87491  17.09606 


INTERCORRELATICN   MATRIX 


COLLECTOR    CURRENT 

1 

1.00 

BETA 

2 

0.29 

1.00 

VBE 

3 

0.86 

0.67 

1.00 

LNUC,    MICROAMPS) 

4 

0.86 

0.68 

1.00 

1.00 

LNIIC,     MICR0AMPS)**2 

5 

0.90 

0.60 

1.00 

0.99 

1.00 
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C3  MICROAMPERE  LINEAR  REGRESSION  PGM,  C3  IC-NOMINAL  DATA,  VCE  =  2.0  VOLT 

CEPENOENT  VARIABLE  IS  NUMBER   -  18  -    BETA 

SAMPLE  SIZE  =   14. C        DEGREES  FREEDOM  =   11.0 
R**2  =  0.918        R  =  0.958         STD  ERR  OF  EST  =    7.705 
ADJ  R**2  =  C.9C3    AOJ  R  =  C.950     ADJ  STD  ERR  OF  EST  =    8.692 


LINEAR  REGRESSION  EQUATION  COEFFICIENTS 

VAR  NC.     NAME                     REG  WT  BTA  WT  WT  SIG    T-RATIO 

0    CONSTANT  TERM           -554. 67056  -20.60101  68.90291    -8.05003 

4  LNIIC,  MICROAMPS)         172.42553  7.22899  19.99124     8.62505 

5  LNIIC,  MICR0AMPS)**2     -11.17180  -6.58700  1.42151    -7.85908 


UEPENDENT  VARIABLE  IS  NUMBER   -  18  -    VBE 

SAMPLE  SIZE  =   14.0        DEGREES  FREEDOM  =   11. 0 
R**2  =  0.999        R  =  l.OCO         STD  ERR  OF  EST  =    0.001 
ACJ  R**2  =  0.999    ADJ  R  =  l.OOC     ADJ  STD  ERR  OF  EST  =    0.001 

LINEAR  REGRESSION  EQUATION  COEFFICIENTS 

VAR  NO.     NAME  REG  WT  BTA  hT  WT  SIG  T-RATIO 

0    CONSTANT  TERM  C. 54023  17.59311  0.00693  77.93624 

4  LNIIC,  MICROAMPS)  C.C2429  0.89277  0.00201  12.07571 

5  LNIIC,  MICR0AMPS)**2  0.00021  0.10742     0.00014  1.45295 

RESULTS  ARE  IN  RAfc  SCORES 
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3   INSTRUCTIONS  FOR  PREPARING  LINEAR  REGRESSION  PROGRAM 
CONTROL  CARDS 

C  1. TITLE  CARD  WITH  FORMAT  =  18A4,  MAY  START  IN  COL  1,  ALL  OTHER 

C  CCNTROL  CARDS  MUST  START  IN  CD  COL  7. 

C 

C  2. MAIN  PARAMETER  CD  -  PROVIDES  I/C  PARAMETERS  TC  PGM. 

C 

C  A.  COL  7-8  NINO  =  12,  NUMBER  OF  INDEPENDENT  INPUT  VARIABLES 

C 

C  B.  COL  10-11  NVAR  =NINO  ♦  NO.  VAR  CREATED  BY  TRANSFORMATIONS 

C 

C  C.  COL  13  ICOR  =  II,  IF  IC0R=1,  THEN  I NTERCORREL AT  ION  MATRIX 

C  APPEARS  IN  THE  OUTPUTt  OTHERWISE  NO  MATRIX  IS  PRINTED. 

C  IF  NVAR  .GE.  15  SET  ICOR  =  0.  THIS  RESTRICTION  DUE  TO  FMT. 

C  BESL"RE  TO  INCLUDE  A  BLANK  CARD  AS  THE  FIRST  INVERT  CARD 

C  WHENEVER  ASKING  FOR  THE  I NTERCORRELATI ON  MATRIX. 

C 

C  D.COL  15  IPRT=CARRIAGE  CONTROL, SETS  PAGE  SPACING  BETWEEN  EONS 

C 

C  0   MEANS  FOR  PRINTER  TO  SKIP  TWO  LINES 

C  -   MEANS  FOR  PRINTER  TO  SKIP  THREE  LINES 

C  1   MEANS  FOR  PRINTER  TC  SKIP   TO  NEW  PAGE 

C  2   MEANS  FOR  PRINTER  TO  SKIP  TO  NEW  HALF  PAGE 

C 

C  fc.  CCL  17  ICHG=1  TO  TRANSFORM  VARIABLES,  FOR  NC  TRANS,  ICHG=0 

C 

C  F.  CCL  19-27   TST  =  F  FORMAT  NC.  .GE.  MAX.  DATA  VAR  ♦  0.5 

C  TST  USED  BY  PGM  TO  SIGNAL  END  OF  DATA. 

C 

C  3.  DATA  VAR  NAME  CDS,  ONE  NAME  PER  CD.  FCRMAT  IS  bA4. 

C 

C  4.  INPLT  CATA  FORMAT  CD,  EX . ( 16X ,6F 10 . 3 »  =  ENTRY  NEEDEO  TO  REAC 

C  6  SETS  OF  INPUT  DATA.  PARENTHESES  ARE  REQUIRED  IN  ENTRY. 

C  IF  A  DIFFERENT  FORMAT  DESIRED  CHECK  INTO  DATA  TITLES  OF 

C  FORMAT  STMT  NO.  615. 

C 

C  5.  CHANGE  CARDS  DEFINE  THE  VARIABLE  TRANSFORMATIONS  FOR  FORMING 

C  LINEAR  REGRESSIONS  EQUATIONS  AS  DEFINED  BY  THE  INVERT  CDS 

C 

C  A.  TRANSFORMATION  CARD  INPUT  CODES  APE... 

C 

C  I   COLS  7-8     IN  FIXED  PCINT  FORMAT  ONLY 

C  J   COLS  10-11   IN  FIXED  POINT  FORMAT  ONLY 

C  K   COLS  13-14   IN  FIXED  PCINT  FORMAT  ONLY 

C  L  OR  DL  COLS  19-27   IN   FLOATING  POINT  FCRMAT  ONLY 
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c 

C  8.   TRANSFORMATIONS  OPERATIONS  COOE  IS.. 

C  01   A(K)=A<J)*»Dt 

C  02   A(Kl=l./( A(J)**OL» 

C  03   Al  K)=A< J)*A(L) 

C  CM       A<K)=AI J)/A(L) 

C  05   A(K)=A( J)-A< i) 

C  06   A(K)=SQRT(At J)  ) 

C  07   AU)=ELOG(A<  J)  ) 

C  08   AIK)=EXPF( A< J)  ) 

C  09   A(K)=A< J»*A< J)  > 

C  10   A(K)=ICDT 

C  11   A(K)=A(J)*DL 

C  12  A(K>=A(Jl*(DL> 

C 

C  C.  AN  EXAMPLE  TRANSFORMATION  SPECIFICATION  CODE  IS... 

C 

C  INDEX        II  JJ  KK  LLLLLL 

C  CARO  COL    7890123*56789 

C  TRAMS  COOE  C4  03  07  06. 

C 

C  THIS  CODE  INSTRUCTS  THE  PGM  THAT  DESIRED  TRANSFORMATION 

C  IS  04   A(K)=A( JI/AILI.  PGM  IS  TO  DIVIDE  VAR  3  BY  VAR  6  AND 

C  STCWE  THE  RESULT  IN  VAR  AC7I. 

C 

C  D.  NOTE   DL  IS  USED  OCCASIONALLY  AS  A  FIXED  POINT  CONSTANT 

C  MULTIPLIER  OR  EXPONENT.  DL  IS  ENTERED  IN  COLS  8-12   IN 

C  PLACE  OF  L..  NOTE  THAT  L  AND  DL  ARE  NEVER  REQUIRED  IN  THE 

C  SAME  TRANSFORMATION! 

C 

C  E.  ENTRY  00  IN  COLS  7-6  MEANS  MC  MCRE  CHANGE  CARDS  PRESENT* 

C 

C  7.  AN  ENC  OF  DATA  CARD  IS  PUT  AT  THE  CONCLUSION  OF  EACH  SET. 

C  THIS  END  CARD  MUST  HAVE  A  VALUE  (IN  THE  FIELD  OF  THE  FIRST 

C  VARIABLE)  THAT  IS  GREATER  THAN  TST  BY  0.1  OR  MORE,   NOTE 

C  VALUE  FOR  TST  WAS  ENTERED  ON  THE  MAIN  PARAMETER  CARD  2. 

C  SEVERAL  SETS  OF  DATA  CAN  BE  CONCATENATED  IN  ONE  CMPTR  RUN. 

C 

C         8.  INVERT  CARDS  DEFINE  MODELS  USING  COLS  7  -  37, WITH  COL  7  FOR 

C  VAR(l). 

C 

C  A.  AN  EXAMPLE  INVERT  CARC  CODE  IS.... 

C  CARO  COL  7  8  9  0  12  3 

C  ENTRY     2    C  0  1  1  0  C 

C 

C  B.  REGRESSICN  MODEL  CODES.... 

C 

C  COL  7=C  =  VARIABLE  NOT  USED 

C  1  =  INDEPENDENT  VARIABLE 

C  2  =  DEPENDENT  VARIABLE 

C  3  =  RESTART  P6M  WITH  ALL  NEW  VARIABLES,  CAUSES 

C  PGM  START  THRU  THE  SET  OF  MODELS  ONCE  AGAIN  AS  DEFINED  BY 

C  THE  INVERT  CD  USING  NEW  DATA. 

C  COL  7=4  =  END  OF  PGM  CARD,  TERMINATES  ALL  PROGRAM  CALC. 

C  3  =  RESTART  PGM  WITH  ALL  NEW  VARIABLES 

C 

c 
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9.       SAMPLE    CATA    DECK 


C] 

.    MICROAMPERE    LINEAR    REGRESSION    PGM, 

CI    TRANSISTCR    OATA,    VCE    =    2.0    VOLTS 

c 

03    05     1    2    1 

12000.0 

NIND.NVAR, 

ICOP.IPRTt 

ICHG.TST 

c 

COLLECTOR    CURRENT 

VARIABLE    NAME    CARD 

c 

BETA 

c 

VBE 

c 

LMIC,     MICROAMPS) 

c 

4.N1IC,    MICR0AMPS)**2 

c 

(  16X.3F 10.3  J 

INPUT    DATA 

FORMAT    CC 

c 

07    01    04    0. 

0 

VARIABLE    TRANSFORMATION    CD 

c 

03    04    05    4. 

C 

c 

00 

c 

3.00 

C.00025C 

38.300 

90.667 

0.632 

2.368 

0.000 

C1DAT    99 

c 

3.00 

Oi.000500 

62. IOC 

95.200 

0.652 

2.348 

0.000 

C1DAT100 

c 

3.00 

O.C01000 

112. 3CC 

100. 40C 

0.672 

2.328 

0.000 

C1DAT101 

c 

3.00 

0.002000 

213. 5CC 

101.200 

0.691 

2.307 

0.002 

C1DAT102 

c 

3. CO 

0.C03000 

320.800 

107. 30C 

0.703 

2.297 

0.000 

C10AT103 

c 

3.00 

O.C04000 

428.300 

107.500 

0.710 

2.290 

0.000 

C1DAT104 

c 

3.00 

0.CO500C 

535.300 

107.000 

0.716 

2.281 

0.003 

CIDAT105 

c 

3.00 

C.C0600C 

642. 4CC 

107.100 

0.720 

2.279 

0.001 

C1DAT106 

c 

3.00 

0.007000 

748. 80C 

106.400 

0.724 

2.274 

0.002 

C1DAT107 

c 

3.00 

C. 008000 

854. ecc 

106. COO 

0.728 

2.272 

C.000 

C1DAT108 

c 

3. CO 

0.C1OOOO 

1C79.000 

112. IOC 

0.734 

2.266 

0.000 

CIDAT109 

c 

3.00 

C.01200C 

1287.999 

104. 50C 

0.737 

2.263 

0.000 

C1DATU0 

c 

3.00 

C.01500C 

1589. OCO 

100.334 

0.746 

2.254 

0.000 

C1DAT111 

c 

3.00 

C.C20000 

2113. OCO 

104. BOO 

0.755 

2.245 

0.000 

C1DAT112 

c 

3.00 

0.025000 

2620.000 

101.400 

0.762 

2.238 

C.000 

C1DAT113 

c 

3.00 

C. 030000 

3117. 999 

99.60C 

0.767 

2.232 

0.001 

C1DATU4 

c 

3.00 

O.C3500C 

3608.000 

98.000 

0.772 

2.226 

0.002 

C1DAT115 

c 

3. CO 

0.040000 

4094.999 

97.40C 

0.776 

2.224 

0.000 

C1DAT116 

c 

3.00 

C.05000C 

5054.996 

96.C0C 

0.782 

2.217 

0.001 

C1DAT117 

c 

3.00 

C.C60000 

5929.996 

87.500 

0.790 

2.210 

0.000 

CIDAT118 

c 

3.00 

O.C70000 

6e64.996 

93.500 

0.794 

2.206 

0.000 

C1DATU9 

c 

3.00 

0.080000 

7666.996 

80.20C 

C.800 

2.200 

C.000 

CIDAT120 

c 

3.00 

0    000100 

24.700 

0.0 

0.603 

2.397 

0    000 

CIDAT    98 

C999999999999999999999999999999999999999TST    CARD    TU    SIGNAL    END    CF    DATA 

C 

COMPLETELY 

BLANK    CD 

SINCE    IC0R= 

1 

C 

C 

2    0    11 

REGRESSION 

MCCEL    OEFN    -INVERT    CD 

c 

0 

0    2    11 

c 

4 

END    CF    PROGRAM    CARD 

c 

c 

****************END    CF    SAMPLE    DATA    DECK********************* 

TJTPICAL     INPIT    OAT*    DECK 
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C3  MILLIAMPERE  LI 
03  05  1  2  1 
COLLECTOR  C 
BETA 
VBE 

LNilC,    HILL 
LNdCt    HILL 
I  16*, 3F  10.3 
07    01    04    0. 
03    C4    05    4. 
00 
2.00       0.C010C0 
2.00      0.005000 
2.00       0.010000 
2.00       0.C12000 
2.00       0.015000 
2.00       0.020000 
2.00      0.025000 
2.00       0.030000 
2.00      0.035000 
2.00      0.C40000 
2.00       0.050000 
2.00       0. 060000 
2.00       0.C7000C 
2.00      0.C80000 


NEAR    REGRESSION    PGM,    C3    TRANSISTCR    DATA,    VCE    *    2.0    VOLTS 

12000.0 
URRENT 


IAHPS1 

IA^PS)**2 

) 

0 

C 


99999999999999999 


C. 106500 
C.5400C0 
1.080000 
1.292999 
1.615999 
2.139999 
2.656000 
3.183000 
3.688000 
A. 191999 
5.212999 
6.176000 
7.207000 
8.230000 
9999999999 


0.0 
108.375 
108.000 
106.500 
107.667 
104.600 
103.600 
105.000 
101.000 
100.800 
102.100 
96.300 
103.100 
102.300 


0.658 
0.702 
0.720 
0.725 
0.730 
0.738 
0.745 
0.749 
0.754 
0.758 
0.763 
0.770 
0.772 
0.775 


999999999999  9999999999999 


.342 
.298 
.280 
.275 
.270 
.262 
.255 
.251 
.246 
.242 
.237 
.230 
.228 
.225 
999999 


0.0 

0.000 

C.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 


C30AT  27 

C3DAT  28 

C30AT  29 

C3DAT  30 

C30AT  31 

C3DAT  32 

C3DAT  33 

C3DAT  34 

C3DAT  35 

C30AT  36 

C3DAT  37 

C30AT  38 

C3DAT  39 

C30AT  40 


999999  9999999999999999 


0  2  0 
0  0  2 
3 


1     1 
1     1 
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C3   X 


2. 
2. 
2. 
2. 
2. 
2. 
2. 
2. 
2. 
2. 
2. 
2. 
2. 
2. 
9999 


ICROAMPERE    L 
03   05    1    2 
COLLECTOR 

BET* 
VBE 

LNUC,  MIC 
LNIiCt  NIC 
I  16K.3F10. 
07  CI  04  0 
03  C4  C5  4 
00 
00  C.C01C0C 
0.005000 
O.C10000 
O.C12000 
0.C1500C 
0.020000 
O.C25000 
O.C30000 
0.035000 
0.04000G 
0.050000 
C. 060000 
0.070000 
C.08000C 


INEAR  REGRESSION  PGM.  C3  TRANSISTOR  DATA*  VCE  *  2.0  VOLTS 

1  12000.0 

CURRENT 


ROAMPS) 

R0AMPS)**2 

31 

.0 
.0 


999999999999 


1C 
54 
108 
129 
161 
213 
265 
318 
368 
419 
521 
617 
720 
822 
99999 


6.500 
C.000 
C.OCC 
2.999 
5.999 
9.999 
8.0CC 
3.0C0 
8.000 
1.996 
2.996 
5.996 
6.996 
9.996 
99999 


0.0 
108.375 
108. OOC 
106.500 
107.667 
104. 800 
103. 60C 
105. OOC 
101.000 

ioo. eoc 

102.100 

96.300 

103. IOC 

102.300 

99999999999 


0.658 
0.702 
0.720 
0.725 
0.730 
0.738 
0.745 
0.749 
0.754 
0.758 
0.763 
C.770 
0.772 
0.775 


1.342 

0.0 

1.298 

0.000 

1.280 

0.000 

1.275 

0.000 

1.270 

0.000 

1.262 

0.000 

1.255 

0.000 

1.251 

0.000 

1.246 

0.000 

1.242 

0.000 

1.237 

0.000 

1.230 

C.000 

1.228 

0.000 

1.225 

C.000 

C3DAT  27 

C3DAT  28 

C3DAT  29 

C30AT  30 

C3DAT  31 

C30AT  32 

C3DAT  33 

C3DAT  34 

C3DAT  35 

C3DAT  36 

C30AT  37 

C3DAT  3B 

C3DAT  39 

C3DAT  40 


9999999999999999999999999999999999999999999 


C  2  0 
0  C  2 
5 


1     1 
1     1 
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5   LINEAR  REGRESSION  INPUT  DATA  TRANSFORMATION  PROGRAM 

THIS  PROGRAM  CALCULATES,  USING  MEASURED  CARD  INPUT  DATA,  BETA, 
MEAN  COLLECTOR  CURRENT, ICM,  DATA  FORMAT  CONVERSION, 
AND  PUNCHES  CARDS  WHICH  ARE  THEN  USED  AS  INPUT  OATA  TO  LINEAR 
REGRESSION  MODEL  COEFICIENT  CALCULATING  PROGRAM. 

COMMON  VCEI26C),I8(260>,IC<260),VCB!2601,VBE!26C) 
COMMON  BETA!260),TVCE(260)  ,-VCEDE  V  <  260  ) 
COMMON  OIC!260),DIB(26C»,ICAI26C),ICM!260I 
200  FORMAT!*  •  ,4X  ,F9.  3,  1X,F10.6,  IX, F  10.  5,  1X,F6.3,5X  ,F6  .  3  ) 

300  FORMAT! • 1»» 

301  FORMAT! ' 2' ,2X , • VCE • ,6X, • IB ' , 5X, • IC-MEAN* ,5X , 'BETA • , 7X, • V8E«  ,  7X , 
1    •VCB',6X,«VCEDE  W/l 

400F0RNATI"  •,F5.2,2F10.6,F10.3,3F10.3,6X,»C10AT«, 131  MA.REGSN 

400F0RMAT!'  ',F5.2,  F 10. 6, F 10. 3, F  10. 3, 3F10. 3,6X, 'C1DAT • ,  13  )  MICR.REG 

401F0RMATM  «,F5.2    , 2F  10. 6, F 10.3 ,3F 10. 3,6X, 'C10AT • , I  3 ,5X , 3F12.3I    MA.REGSN 

401  FORMAT! •  • ,F5 .2 , 1 F 10 .6,F 10 . 3  ,F 10. 3, 3F 10. 3 ,6X, 'C ID AT • , I  3, 5X.3F 12 .3  IMICR. REG 

402  FORMAT! •-•, 22 X • **********  ENO  OF  DATA  *•*********•) 
REAL  IB.ICICA,  ICM 

N=000 

BETAI1)    =    C.CCC 
DIB! 11=0. 
Old  I)    =    0. 
VBE!  I)    =    O.C 

icm  =  ccoocooi 

IB(1)  =  CCCCCOCl 
VCE«1)  =  C.5 
WRITE(6,3CC) 
DO  a3  J=2,261 

12  N=N*1 

BETA(J)  =  C.000 
10  READ!5,20C,ENC=14)VCE! Jl ,  IB« J  )  , ICIJ I , VBE IJ> ,VCB!J) 

13  CONTINUE 

14  00  16  J=2,  N 

IF(VCE!J«1)  .NE.  VCEIJI)  GOTO  15 
8  DIB!J*1)  =  IBIJ+1)  -  IB! J } 
DIC4J  +  1)  =  IC(  J  +  l)  -  ICIJ) 
IF!IB!J*ll  .EG.  IB! J  )  )  GOTO  15 
ICM!J*ll  =  !IC!J*1)  ♦  IC!J))/2.C 
ICM!J*1)  =   ICM! J»l)*1000.0 
BETAIJ+1)  =  DIC(  J  +  11/DIBIJ-U) 

15  TVC£( J)=VBEi J)  ♦  VCB1J) 
VCEOEVIJI  =  VCEIJI  -  TVCEiJI 
ICA4J)  ■  (  ICM)  1/1.0 

16  CONTINUE 
VCE«N*1I  =  1CC.C 

30  DO  40   J=2,  N 

K=J-1 

IFIVCEiJ+1)  .GT.  VCEU))  WRITE!6,301I 

WRITE16,401)VCE( Jl, IB!J),ICN! J), BETA! J) ,  VBEIJ) , VCB! J  )  , VCEDEVi J  I , K 
l.BETAIJ   J.CICIJ   J.DIBIJ   I 

WRITE  I  7,400 )VCE< J), IB! J), I  CM  I  J) .BETA! J) , VBE I  J) , VCB( J  I ,VCEDEV! JI,K 
40  CONTINUE 

V.RITEI6.4C2) 

STOP 

END 
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APPENDIX  0 

REGRESSION  MODEL  EFFECTIVENESS  PROGRAM 

1   REGRESSION  MODEL  EFFECTIVENESS  PROGRAM  LISTING 

COMMON   VCEI200I,  BET AMC20CI , VBEM1200 I , BETAC 1200 ) ,VBECI 200  I 

COMMON   VDIFF(200I,BDIFF!200>,B<5),V!5),K,M,TITLEI8I .UNITS! 5) 

COMMON   INT 

REAL   ICMI2CC1 
102  F0RMAT(6X,8A4) 
104  F0RMAT(6X,4F12.5,21X,5A1) 
106  FORMAT!6X,4F12.5,21X,5Al) 
108  FORMAT! I  1  ,F5  .2  .  10X , 3F 10. 3» 

200  FORMAT! M'i  5X, 'ANALYSIS  OF  REGRESSION  EQUATION  FOR  BETA  -  «8A4/ 
l«-«  ,5X,«BETA  =  ',F10.5,«  ♦  '.F10.5,'  *  LN  <IC)  ♦  ',F10.5,«  *  • 
2*  !LN!IC) )**2',/   ) 
202  FORMAT! '0«,7X,'VCE'  .  1 IX, ■ IC« . 10X, • MEASURED*  ,6X, 'CALCULATED* 8X, 

1  'BETA', /,36X,«BETA«,11X, 'BETA', 8X, 'MEAS-CALC  ,/, 

2  7X,'V0LTS',6X,  '  (  •  , 5AI , • AMPS ) • , 7X , • !  —  )  •»11X,«C  —  ) '  ,  11X,' 1  —  )•     I 

204  F0RMATI6X,F6.2,4F15.5) 

205  F0RMATM1',  5X, "ANALYSIS  OF  REGRESSION  EQUATION  FOR  VBE  -  ',8A4,/ 
l'-'tSX.'VBE  =  ',FS.5,'  ♦  •,F9.5.«  *  LN  1IC)  ♦  ',   F9.5,    •  *  • 

2*  (LNIICn**2(t/  ) 

206  FORMAT! 'O'.TX,  'VCE'  , 1 1 X, • IC ', 10X, • MEASURED •  ,6X, 'C ALCULA TED' 8X , 

1  'VBE* ,/,37>, 'VBE  '.10X,'  VBE  •  ,8X , «ME AS-CALC , /  , 

2  7X,  'VOLTS'  ,6X,  •  !«,5Al,  'AMP  SI  '  ,7X,«  ! — > • , 1 IX , • I -- ) • ,  1 1 X , • I  —  ) •     ) 
208  F0RMAT(6X,F6.2,4F15.5) 

10  K=0 

M=l 

REA0I5,  102,END=40)    ITIILEIII),  11  =  1,8) 

REAC<5,104)  Bi  1  ),B(2),B(3) ,B<4)  ,IUNITSIKK),  KK=1,5) 

REA0(5,106)  VI 1), VI2).V! 3) ,V!4)  .1UNITS1KK),  KK=1,5) 

BI4)  =  C.C 

VIA)  =  O.C 

DO  15  1=1, 2CC 

REAO15,108,Ef\D=2C)  IN  T  ,  VCE  (  I)  ,  I  CM  (  I  I  ,  BETAM!  I  )  ,  VBEM!  I  ) 

IFIINT  .GE.  M  GCTC  20 

K  =  K  ♦  1 
15   CONTINUE 
20  DO  30  J=1,K 

BETACIJ)  =  B!l)   ♦  B(2)*(ALUGi  1CM1J))  )  ♦  B  I  3  )  *  (  <  ALOG  <  I  CM  (  J  )))  **2  ) 
VBECIJ)   =  V!l)  ♦  V!2)*!AL0G!ICM(JI)  I  ♦  V! 3) *( ( ALOG < ICM! J ) ) ) **2) 

BDIFF(J)  =  BETAMIJ)  -  EETACIJ) 

VOIFF(J)  =  VBEM(J)  -  VBECIJ) 
3C  CONTINUE 

HRITE!6,200)  ! T  I TLE I  J J  ) , J J=  1 , 8 )  ,  BI 1 )  ,8 1 2  ) , B I  3 ) 

WRITEI6.2C2)   ILNITSiKK),  KK=1,5) 

WRITE(6,204)  ( VCE (  i ) , ICM ( I ) , BET AM  I  I  )  ,  BE TAC (  I  I , BO  I FF i  I) ,  1  =  1, K) 

MRITE(6,205)  I T 1 TLE i J J ) , J J= 1 , 8  I  ,  V (  1 ) , V( 2  )  , V( 3 ) 

WRITEI6,2C6)   iUNITS(KK),  KK=1,5) 

WRITE!6,2C8)  ( VCE! I  )  , ICM! I ) , VBE f ! I » . VBEC I  I )  ,  VOIFF ! I ) ,  1  =  1, K) 

INT  =  0 

GOTO  IC 
40  STOP 

END 
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ANALYSIS  OF  REGRESSION  EQUATION  FOR  BETA  -  C3  IC-NOHINAL  VCE-2.0  VOC 


BETA 


103.31700  ♦ 


18.08138  *  LN  (Id  ♦   -11.17180  *  (LN(ICI)**2 


VCE 

VOLTS 

0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
CO 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 


IC 

MEASURED 

CALCULATED 

BETA 

BETA 

(MILL1AMPS) 

(  — ) 

<--) 

C.1C6C0 

0.0 

6.46486 

0.540CO 

108.37500 

87.93373 

1.08000 

1C8. 00000 

104.64238 

1.29300 

106.50000 

107.22557 

1.61600 

107.66699 

109.42172 

2.14000 

104.79999 

110.60689 

2.66800 

103.59999 

110.31653 

3.18300 

1C5.00000 

109.27560 

3.68800 

101.00000 

107.88643 

4.19200 

IOC. 79999 

106.28398 

5.21300 

102.C9999 

102.71432 

6.  17600 

96.29999 

99.20450 

7.20700 

103.09999 

95.44937 

8.23000 

102.29999 

91.79503 

BETA 
MEAS-CALC 
I  —  J 

-6.46486 
20.44127 

3.35762 
-0.72557 
-1.75473 
-5.80690 
-6.71654 
-4.27560 
-6.88643 
-5.48399 
-0.61433 
-2.90451 

7.65062 
10.50496 


ANALYSIS  OF  REGRESSION  EQl;AT  ICN  FOR  VBE 


C3  IC-NCMINAL  VCE=2.0  VOC 


VBE  = 


VCE 

VCLTS 

0.0 
0.0 
0.0 
0.0 
0.0 
0-0 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 


0.71791  ♦    0.02716 

*  LN  I IC)  ♦ 

0.00021  * 

IC 

MEASURED 

CALCULATED 

VBE 

VBE 

(HILLIANPS) 

I—) 

(  —  1 

0.10600 

C. 65800 

0.65801 

C.  54000 

C. 70200 

0.70125 

1.C8000 

0.72000 

0.72000 

1.293CC 

0.72500 

0.72490 

1.616C0 

C. 73000 

0.73099 

2.140C0 

C. 73800 

0.73869 

2.65800 

0.74500 

0.74466 

3.  183C0 

0.74  900 

0.74964 

3.688C0 

0.75400 

0.75371 

4.19200 

C.7580C 

0.75727 

5.21300 

C. 76300 

0.76333 

6.  176C0 

C. 77000 

0. 76806 

7.207C0 

C. 77200 

0.77237 

8.23000 

C. 77500 

0.77609 

VBE 
MEAS-CALC 

-0.00001 
0.00075 

-0.00000 
0.00010 

-0.00099 

-0.00069 
0.00034 

-0.00064 
0.00029 
0.00073 

-0.00033 
0.00194 

-0.00037 

-0.00109 
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ANALYSIS    OF    REGRESSION    EQUATION    FOR    BETA   -   C3    IC-NOMINAL    VCE«2.0    VOC 
BETA    =   -554.67C41    ♦       172.42562    *  LN    (JO    ♦      -11.17180    *    <LN(IC))**2 


VCE 

IC 

MEASURED 

CALCULATED 

BETA 

BETA 

BETA 

MEAS-CALC 

VOLTS 

IMICROAMPSI 

1—1 

1— 1 

1  — 1 

0.0 

106.50000 

0.0 

6.78583 

-6.78583 

0.0 

540.00000 

108.37500 

87.93408 

20.44092 

0.0 

1C80. 00000 

108. 00000 

104.64282 

3.35718 

0.0 

1292.99878 

106.50000 

107.22583 

-0.72583 

0.0 

1615.99878 

107.66699 

109.42187 

-1.75488 

0.0 

2139.99878 

104. 79999 

110.60718 

-5.80719 

0.0 

2658.00000 

103.59999 

110.31689 

-6.71690 

0.0 

3183.000C0 

105.00000 

109.27612 

-4.27612 

0.0 

3688.CC0C0 

101.00000 

107.88672 

-6.88672 

0.0 

4191.  99219 

IOC. 79999 

106.28442 

-5.48444 

0.0 

5212.99219 

102.09999 

102.71484 

-0.61485 

0.0 

6175.99219 

96.29999 

99.20483 

-2.90485 

0.0 

7206.99219 

103.09999 

95.44971 

7.65028 

0.0 

8229.99219 

.102.29999 

91.79541 

10.50458 

ANALYSIS    OF    REGRESSICN    EQUATION    FOR    VBE    -    C3    IC-NOMINAL    VCE=2.0    VOC 


E  = 

G.54C23  ♦    C.C2429 

*  LN  (  IC)  ♦ 

C.CC021  * 

[LN<  IC) )**2 

VCE 

IC 

MEASLRED 

CALCULATED 

VBE 

VBE 

VBE 

MEAS-CALC 

VOLTS 

(MICRCAMPS) 

<--) 

1  —  1 

l~l 

0.0 

1C6.500C0 

C. 65800 

0.65820 

-0.00020 

0.0 

540.COO0O 

0.70200 

0.70136 

0.00064 

0.0 

iceo.cooco 

C. 72000 

0.72013 

-0.00013 

0.0 

1292.99878 

C. 72500 

0. 72504 

-0.00004 

0.0 

1615.99878 

0.73000 

0.73114 

-0.00114 

CO 

2139.99878 

C. 73800 

0.73885 

-0.00085 

0.0 

2658.00000 

C. 74500 

0.74482 

0.00018 

0.0 

3183.CC0C0 

0.74900 

0.74980 

-0.00080 

0.0 

3668.00000 

C. 75400 

0.75388 

0.00012 

0.0 

4191.99219 

C. 75800 

0.75744 

0.00056 

0.0 

5212.99219 

C. 76300 

0.76351 

-0.00051 

0.0 

6175.99219 

0.7700C 

0.76824 

0.00176 

0.0 

7206.99219 

C. 77200 

0.77256 

-0.00056 

0.0 

8229.99219 

C. 77500 

0.77629 

-0.00129 
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PROGRAM  COMMENTS 


1.  INPUT  CARDS 

A.  TITLE  CARO  -   READ  IN  TRANSISTOR  TYPE,  METHOD  OF  OBTAINING 
IC  -  WHETHER  NOMINAL  VALUE  OR  MEAN  VALUE  OF  IC  IS  USEO  IN 
THE  COLLECTOR  CURRENT*  IC,  DATA,  ANC  THE   COLLECTOR  TO 
EMITTER  VOLTAGE,  VCE. 

8.  BETA  EOLATION  CARO  -   COEFFICIENTS  CF  THE  REGRESSION 

MODEL  BETA  EQUATION  ARE  READ  IN  BY  FORMAT (6X.4F 12.5, 5A1) . 
NUMERIC  VALUES  ARE  OBTAINED  FROM  REGRESSION  PROGRAM  OUTPUT. 
THE  LAST  5A1  CHARACTERS  CCMPRISE   THE  COLLECTOR  CURRENT 
UNIT'S  PREFIX  OR  MULTIPLYING  FACTOR.  THAT  IS  MICRO  IMPLIES 
MICROAMPERE  UNITS,  MILLI  IMPLIES  MI  LI  I  AMPERE S  UNITS,  AND 
5  BLANKS  SIMPLY  MEAM  THAT   THE  CURRENT  UNITS  ARE  IN 
AMPERES.  THE  PREFIX  READ  IN  IS  PLACED  IN  THE  OUTPUT  UNDER 
THE  IC  HEADING. 

C.  VBE  EQUATION  CARD  -   COEFFICIENTS  OF  THE  REGRESSION 
MODEL  VBE  EQUATION  ARE  READ  IN  BY  FORMAT (6X ,4F 12. 5, 5A1 I . 
NUMERIC  VALUES  ARE  OBTAINED  FROM  REGRESSION  PROGRAM  OUTPUT. 

D.  OATA  CARDS  -  READ  USIKG  FORMAT ( 1 1 ,F 5.2 ,3F 10 .3 ) .  CARDS 
MAY  BE  OBTAINED  FROM  THE  LINEAR  REGRESSION'S  PROGRAM'S 
PUNCHED  OUTPUT  DATA  SET. 

E.  END  OF  DATA  CARD  -  ANY  INTEGER  GREATER  THAN  I  ENABLES 
PROGRAM  TO  BEGIN  READING  ANCTHER  DATA  SET.  IF  A  /*  CARD  IS 
EMCOLNTERED,  THEN  CONTROL  IS  RETURNED  TO  SYSTEM  MONITOR. 
INTEGER  MUST  BE  PLACED  IN  COLIMN  1.  PROGRAM  WILL  NOT  WORK 
WITHCLT  THIS  CARD. 


2.  SAMPLE  CATA  DECK 

CI  IC-NOMINAL  VCE=2.C  VDC  MILLIAMPERES 

CI  2.  B  MILLI 
CI  2.  V  MILLI 


104.48077 

-5.28263 

-3.90C85 

0.74633 

C. 02786 

-0.C0499 

0.007 

O.C 

0.404 

C.017 

99.50C 

0.604 

C.029 

84.20C 

0.631 

C.052 

92.120 

0.651 

0.101 

96.880 

0.669 

0.198 

97.70C 

0.688 

0.300 

101. 50C 

0.699 

0.403 

102.900 

0.708 

C.808 

101.700 

0.727 

1.012 

1C2.000 

0.734 

6.538 

90. IOC 

0.795 

7.293 

75.5CC 

0.802 

3   IA  NO.  .GR.  1  IN  COL  1  PERMITS  RUN  OF  ANOTHER  DATA  SET  OR  TERMINATES  PGM  RUN) 
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CI    IC-NOMINAL    VCE=2.0    VDC  MICROAMPERES 

CI    2.    B   MICRO 
CI    2.    V   MICRO 


-45.16509 

46.60955 

-3.90085 

0.31567 

0.09663 

-0.C0499 

6.750 

0.0 

0.404 

16.700 

99.50C 

0.604 

29.330 

84.200 

0.631 

52.360 

92.120 

0.651 

706.300 

102.000 

0.724 

eC8.C0C 

1C1.70C 

0.727 

1011.999 

1C2.0CC 

0.734 

5636.996 

83.50C 

0.792 

6537.996 

90.100 

C.795 

7292.996 

75.50C 

0.802 

END    OF    SAMPLE    DATA    DECK 


T.VPICAL     INPLT    DATA    CECK 
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C3  IC-NOMINAL  VCE=2.0 

VDC 

103.31701 

18. 06139 

-11.17180 

0.71791 

0.02716 

0.08021 

C.106 

O.C 

0.658 

0.540 

108.175 

0.702 

1.080 

108.000 

C.720 

1.293 

106. 50C 

0.725 

1.616 

107.667 

0.730 

2.140 

1C4.800 

0.738 

2.65e 

1C3.60C 

0.745 

3.183 

1C5.00C 

0.749 

3.686 

101. OOC 

0.754 

4.192 

100.800 

C.758 

5.213 

102.100 

0.763 

6.176 

96.900 

C.770 

7.207 

103.100 

0.772 

8.230 

1C2.30C 

0.775 

C3  IC-NOMINAL  VCE=2.0 

VOC 

-554.67C56 

172.42553 

-11.1718 

C. 54023 

C.  02429 

0.00021 

1C6.50C 

O.C 

0.658 

540. COO 

108.375 

0.702 

1C80.000 

108. COC 

0.720 

1292.999 

106. 50C 

0.725 

1615.999 

107.667 

0.730 

2139.999 

104.800 

0.738 

2658.000 

103. «0C 

0.  745 

3183. OOC 

105. OOC 

0.749 

3688. OOC 

101. OOC 

0.754 

4191.996 

ico.eoc 

0.758 

5212.996 

102. ICC 

0.763 

6175.996 

96.30C 

C.770 

7206.996 

103. IOC 

0.772 

8229.996 

102.300 

0.775 

MILLIAMPERES 
C3  2.  B  MILLI 
C3  2.  V  MILLI 


MICROAMPERES 
C3  2.  B  MICRO 
C3  2.  V  MICRO 
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NLAP  OUTPUT  SUBROUTINE  LISTING 
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12 
11 

14 


SUBROUTINE    ECB25 

COMNON    CCSAV(6,2C0) ,  IUWORDI 74) 
DOUBLE    PRECISION    X(200) 

COMMON   NMAX, NNODE, N TERMS, NUMBL,NUNBR,NUM8C, I RTN.NTR AC E.NSWTCH.KTO 
1    NPRINT(IO) 
COMMON    E(2CC),EMIN1200),EMAK(200) ,  AMP ( 200  ) , AMPM IN < 200 » , AMPMAX < 200 
COMNON    Y(200) , YM INI  200 ) , YMAX< 20C) ,NINIT(200) ,NF IN ( 200 ) , MODE  1 ( 200 ) 
COMMON    Y TERM « 200 ) , YTERMHI 200 ) , YTERMLl 200) , I  ROW T 1 200) ,  IC0LTI200) 
COMMON    ERROR1,ISEQ,MSEQ,MO,«UMNO,VFIRST<50) ,  VSECNOI 50 ) , VL AST ( 501 
COMiON    MOB  AN ( 5C) ,  MOP  ARM ( 50 ) , MCS TE P < 50), IwCCUT(A) 

THE    FOLLOWING    VARIABLES    ARE    USED    ONLY    IN    THE    ECAP    D.C.    ANALYSIS 

COMMON    AX1.SMLEPI50) ,CURR ( 200 ) , SMLE ( 20C ) ,EQlXUR I  50 ) ,EX(200) 
COMMON    EB(200),AMPXI200),AMPB<2  00) , VNOMI 50) , STDSQI 50 ) ,L  ,M,ITOL 
COMNON    JX1,JX4,  JX5, DELTA  ,DUM1<  28) 

COMMON    LANGU217)  ,NTR,LANG2<47) 

COMMON    MATAI2C0,4,3),YX(200), YB(200) , YTERMX ( 200 ) , YTERMB ( 200 ) 
COMNON    HCMAX(50)tUCMIN(50) 
COMMON       CCSAV(2,200) 

OOUflLE    PRECISION    SMLEPiCURR , SMLE , EQOCUR 

IF(NTRACEI3,4,3 
F0RNATI6H    ECB25    ) 
WRITEI6,21 

OUTPUT    NODE    VOLTAGES 

KNNODE    =    NNODE    ♦     1 
DO    5    1=1, NNODE 
CCSAV<6,I)=SMLEP< I ) 
DO    6    K=KNNODE.    NMAX 
CCSAVI6.KI     =       .0 

BRANCH    VOLTAGES 

00   10     1=1, NMAX 

SMLE(I)    =    0.0 

J-MNITI  I  I 

IMJ  111,  11,12 

SMLE I  I)    =    SMLEPCJ) 

K=N* INJ I) 

IF    IK    I     1C,     10,     14 

SMLEIII  =  SMLEII)  -  SMLEP(K) 


0C250000 

0C250010 
,DC250020 
DC250030 
JOC250040 
DC250050 
DC250060 
DC250070 
DC250080 
DC250090 
DC250100 
OC25OH0 
0C250120 
DC250130 
DC250140 
DC250150 

DC250170 
DC250180 
DC250190 

DC250200 
DC250210 
DC250220 
DC250230 

DC250250 
DC250270 
DC250280 
DC250290 
DC250292 


DC250293 
DC250294 
DC250400 
OC250410 
DC250420 
DC250430 
DC250440 
DC250450 
DC250460 
DC250470 
DC250480 
DC250490 
DC250500 


10    CONJINUE 

00    16    IM.NMAX 

16  CCSAVU,  I  )«6*LE<  I) 

ELEMENT    VOLTAGES 

DO    17    IM.NPAX 

SML£ ( I)»SMLE(  I  )*EX{  I  ) 

17  CCSAVI 3, I  )>SMLE< I) 

ELEMENT    CURRENTS 

103    00    20     IM.NHAX 
20   £UR*( I)=YX(l )*SfLE( I ) 
IF(NTERMS)21. 22, 21 


21 


24 
23 
22 
21 


DO    23    I-1,NTERMS 
NR*IROWT(II 

NC=ICOLT<I ) 

IF<SMLE(NC))2<.,  23,2<t 

CURR(NR)=CURR(NR)*YTERMX< I >*SNLE(NC) 

CONTINUE 

DO    27    I=1,\HAX 

CCSAVU.I  )  =  CURR<  I  ) 


BETA  SAVE 


28 


DO  28   N  =  l.NMAX 
NFBR  =  ICOLT(N) 
NTBR  ■  IROWT(N) 
CDSAVI  1  .N  )  =  0.0 
CDSAV(2,N)  ■ 
CDSAVI 1,N)  = 
IFHCOLT(N) 
CDSAV(2,N)  = 
CONTINUE 


33 


35 
34 
37 


0.0 

YTERM(N) 
EQ.  0)  GOTO  26 
(CDSAVI  1,N) )/( Y(NFBR)) 


BRANCH  POWER  LOSSES 


DO  29  I=1,NMAX 

29  CCSAV(5,I >=CLRR(I>*SMLE<I) 

BRANCH  CURRENTS 

105  DO  30  I=1,N*AX 

30  CURRd  J=CURR(  I)-AMPX«  I  ) 

CHECK  UNBALANCES 


DO  33  I=1,NNGDE 

XI  I)  =  C. 

00  36  I=1«NHAX 

J  =  MNIT<  I) 

K=NFIN( I  I 

IFIK)34,34,35 

X(K)  =  X(K  )*CLRR(  I ) 

IFU)36,3£,37 

X<  J)  =  X(J)-Ct£RR<  I  I 


116 
0C250510 


DC250630 
DC250640 
DC250650 


DC250780 
DC250790 
DC250800 
DC250810 
DC250820 
0C250830 

DC250840 
DC250850 
0C250860 
DC250870 
0C250880 
DC250890 


DCB25092 
DCB25093 
0CB25094 
DCB25095 
DCB25096 


DC251010 
DC251020 
DC251030 


DC251140 
DC251150 
DC25U60 
DC251170 
DC251180 
DC251 190 
DC251200 
DC251210 
DC251220 
0C251230 
DC251240 
DC251250 
DC251260 
DC251270 
OC251280 
0C251290 
DC251300 


117 


36 

38 

40 

141 

142 


106 
42 

398 
399 


400 


401 

402 
403 


404 


405 


CONTINUE 

SUM*0. 

£0   38    I>1*NW30E 

SUM   -    SUN    ♦    OABSI    X(     t    II 

IF ISUM-ERR0R1 1 106 , 106,40 

WRITEI6,141) 

WRITE    16,142) 

FORMAT!    //    43H    SOLUTION   NOT   OBTAINED    TO   OESIREO    TOLERANCE//! 

FORMAT!    6*i    NCCE  S  ,  15X  ,  19H    CURRENT    UNBALANCES    /) 

KMAX-NNODE 

IND-6 

GO   TO    100 

00    42    I»1,NNAX 

CCSAVI2,II*CURR(I) 

IF    4JX4    .GT.    C)    GOTO    900 

GO    TO    1399, 399,399, 900), NTR 

CALL    BCDBN1I3,1,IDW0RDI4)     ) 

IF(IOWORO(5)    .GE.    1)    I0HOR0(73)    ■    0 

IDWORDI731  ■  I0W0RDI73I  ♦  1 

IDWORDI73I  ■  1IDWOROI731  ♦  l)/2 

WRITE (6, 8011  I IOWOR0IIK),  IK"7, 72 ) , IOWORO ( 3) , I  DWORD! 4 ) , I  DWORD  J  73 ) 

WRITE  (6,8001 

IF  (NMAX-NNGCE)  4C0,5C0,401 

MAX«NNODE 

MIN*NMAX 

GO  TO  402 

MAX=NMAX 

MIN=NNODE 

00  403  I=1,KMAX 

XI  I |= 1 .0/Y<  N 

WRITE  (6,e5C)  I  I,E(  I), AMP! ll.XI 1 1 , ICCSAVIJ , I ) ,J  =  l,6) ,1,1  =  1, MINI 

IF  JMIN-MAX)  4C4,«5CC,900 

MIK=MIN*1 

WRITE  (6,8511!  I, El  I), AMP! I ),X(I  ),ICCSAV( J, I),  J  =  l,5),  1=  MIN.MAX) 

IF  IIDW0RCI74)  .LE.  0)  GOTO  900  TEST 


DC251310 
DC251320 
DC251330 
0C251340 
DC251350 
DC251360 
DC251370 
OC251380 
DC251390 
0C251400 
DC251410 
DC251420 


DCA25104 
DCA25105 


DCA25108 


DCB25U7 


DCB251 

DCA25121 

DCA25122 


WRITE  (8,852)  (lOWORO(I),  1-1,  73) 

WRITE!  7,652)  (ICWORD(I),  1*1,  73) 

WRITE  (9,8521  (IDWORDII),  1*1,  73) 

WRITE! 10,852)  (IOWORDII),  1*1,  73) 

WRITE (  1  1,8521  (IDWORDII),  I>1,  73) 

WRITE  I  7, 853*  ( I  DWORD  I  3) , ID WORD! 73) , I ,E (  I ) , AMP! I ) ,x( I  I  ,  ( CCSAVI J, I ) , 
1  J=l,6),  1=1,  MAX) 

WRITE (8, 8  54)  ( I  DWORD ( 3 ) , I DWORCI 73  )  ,  I  ,E 1  1 1, AMP ( I ) , X  I  I  )  , ICCSAVIJ, I ) , 
1  J=i,6),  1=1,  MAX) 

WRITE  1 9, 8 55)  ( I DWCRD ( 3 ) , IOWCRC! 73 ) , I ,E(  I), AMP (I ) , X  I  I ) , COS AVI  1 , 1 ) , 
1CCSAVI  I,  1  1, CCSAVI 2,  I  >, CCSAVI 5,1  ),CCSAV(6,  I) ,  1  =  1,  MAX  ) 

WRITE  I  10,8  551  (  1CW0R0I 3) , IOWOROI 73) ,I,EII),AMP( I ) ,X ( I ) , CDS AVI  2, 1 ) , 
1CCSAVI 1, II, CCSAVI2, II, CCSAVI 5,1 1 .CCSAVI6, I )  ,  1  =  1,  MAX  ) 

WRITE  1 11, 855 1  I  I  DWORD! 3  I , I  DWORD ( 731  ,  I  ,E({),AMPI I ) , X  I  I ) , CDSAVI 2 , 1 ) , 
1CCSAVI 1,1  ),CCSAV(2,I) ,CCSAVI5,I  ),CCSAVI6,I)  ,  1=1,  MAX  I 

GO  TO  900 

800  FORMAT  !•  BRANCH ■ ,3X, 'VOL  I  AGE ', 5X, • CURRENT ' ,7X ,' BRANCH' , 7X ,  DCB25121 
1*  BRANCH  •,  8  X,'eRANCH',7X,'  ELEMENT  ',7X,  'ELEMENT ',  5X  ,•  EL  EMENT  •,  7X  DCB2512? 
2'N0DE',6X,'NODE'/'  NUMBER •, 3X, • SOURCE « ,6X , 'SOURCE •, 6X ,• RES  I  STANCE 'DCB25 12  3 
3, 5X,« VOLT  AGE ',7X, 'CURRENT', 6X,' VOLT AGE ',7X,' CURRENT «,4X,  DCB25124 
4 • POWER  LOSS », 4X, • VOLT AGE «,3X,' NUMBER '/10X,' I  VOLTS) • ,6X , • 1  MA) • ,9X,  0CB25125 
5MOHMS)',7X,  '1V0LTS)',8X,'1«IA)',8X,  •  I  VOLTS)  ',8X,'(MA)  •  ,7X,  DCB25126 
6' I  WATTS  )',6X,' I  VOLTS)  •//,/)  DCB25127 

801  FORMAT! • 1' ,4OX,66Al/'0',40X, 'PRCBLEM  NO.  = • ,2X , A  I , A  1 , 10X, ' SOLUT ION 


118 

1  NO.  «  '  13,////) 
650  FORMAT!'  ' ,4 4.2X.F 10 .3 t    3PF  1  2  .4,  3X  ,  0PF12.2  ,  2X  ,F  1  1 .6  ,  2X  ,  3PF  1  2  .6, 
12X,CPFll.fc,PX,3PFl2.6,2X,0PF 10.4 , 2X ,OPF U .6 , 2X , I  5 ) 

851  FORMATI •  '  ,  14,  2X  ,  F  10. 3  ,    3PF 12 .4, 3X,0PF12 .2 ,2X ,F 1 1 .6 , 2X , 3PF1 2 .6 , 
12X,OPF11.6,2X,3PF12.6,2X,0PF10.4) 

852  FORMAT  (  1  X  ,  •  C  '  ,  1  X , 72A I , I  3 » 

853  F0RMATI'P',Al,«-S',I3,«-B»,I2t    F7  .3  ,3PF8. 2.0PF 1 1 . 2,F7. 3, 3PF8.3, 
10PF7.3,3PF8.3,0PF6.3,0PF7.3I 

854  FORMATI'  P •  ,  Ai , '-S« , I  3 , «-B • , I  2,    F7. 3 , 3PF6.2 ,0PF 11 .2 ,F7. 3, 3PF8.3 , 
10PF7.3,3PF9.3,0PF6.3,0PF7.3I 

855  FORMAT< •  P  •  ,  A  I, «-S • , I  3, «-B • , I  2,    F8. 3, 3PF8.2 ,OPF 1 1 . 2,  F9.3,F8.3, 
13PF9.3,  0PF7.3,CPF8.3) 

856  FORMATI •  P* ,A1 , '-S* , I  3, «-fl • *  I  2, IX , F 7. 3, 3PF8. 2 , OPF 1 1 .2 ,F7.3 , 3PF 10. 3 
1   ,0PF6.3,0PF7.3) 

860  FORMATI  •  •  ,  5X  ,  • P • , A  1 , «-S • , I  3 , '-B'  ,  1 2  ,  5X  ,  •  ECB25  •  ,10X,«  ICCLTP,  13,'  ) 

1  *•»  I4.10X,  MPOhTI  •  ,13, •  )  =«,I4  I 

C  DC251520 

C      OUTPUT  ROUTINE  DC251530 

C  0C251540 

100  LASI=0  DC251550 

150  K*LiAST*l  DC251560 

LAST=LAST*4  0C251570 

IF  ILAST-KMAX)  200,200, 201  DC251580 

201  LAST=KMAX  DC251590 

200  WRITEI6, 2C3)K, LAST, ( X|JI  ,J*K, LAST)  0C251600 

203  F0RMATI1X,I3,1H-, I3.2X.4I 3X,E15.8) )  DC251610 

IF(KMAX-LAST)5CC,5CC,150  DC251620 

C  DC251630 

500  GOTO  106 

900  RETURN  DC251650 

END  0C251660 
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